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Abstract 

Kinetic relations are required in order to characterize nonclassical un- 
dercompressive shock waves and formulate a well-posed initial value prob- 
lem for nonlinear hyperbolic systems of conservation laws. Such nonclas- 
sical waves arise in weak solutions of a large variety of physical mod- 
els: phase transitions, thin liquid films, magnetohydrodynamics, Camassa- 
Holm model, martensite-austenite materials, semi-conductors, combustion 
theory, etc. This review presents the research done in the last fifteen years 
which led the development of the theory of kinetic relations for undercom- 
pressive shocks and has now covered many physical, mathematical, and 
numerical issues. The main difficulty overcome here in our analysis of 
nonclassical entropy solutions comes from their lack of monotonicity with 
respect to initial data. 

First, a nonclassical Riemann solver is determined by imposing a single 
entropy inequality, a kinetic relation and, if necessary, a nucleation crite- 
rion. To determine the kinetic function, the hyperbolic system of equations 
is augmented with diffusion and dispersion terms, accounting for small- 
scale physical effects such as the viscosity, capillarity, or heat conduction of 
the material under consideration. Investigating the existence and proper- 
ties of traveling wave solutions allows one to establish the existence, as well 
as qualitative properties, of the kinetic function. To tackle the initial value 
problem, a Glimm-type scheme based on the nonclassical Riemann solver is 
introduced, together with generalized total variation and interaction func- 
tionals which are adapted to nonclassical shocks. Next, the strong conver- 
gence of the vanishing diffusion-dispersion approximations for the initial 
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value probem is established via weak convergence techniques. Finally, 
the numerical approximation of nonclassical shocks relies on schemes with 
controled dissipation, built from high-order, entropy conservative, finite 
difference approximations and an analysis of their equivalent equations. 
Undercompressive shocks of hyperbolic conservation laws turn out to ex- 
hibit features that are very similar to shocks of nonconservative hyperbolic 
systems, who were investigated earlier by the author. 
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1 Introduction 

1.1 Augmented systems of conservation laws 

Nonlinear hyperbolic systems of conservation laws have the general form 

ii t + f(u) x = 0, u = u(t,x), (1.1) 

with unknown u : R + X ]R — » It, where the given flux / : U — > ]R N is defined on 
a (possibly non-connected) open set U C R N and satisfies the following strict 
hyperbolicity condition: for every v e U, the matrix A{v) := Df(v) admits N 
real and distinct eigenvalues 

Ai{v) < ... < A N (v) 

and a basis of right-eigenvectors ri (c), . . . , Tk/(v). It is well-known that singulari- 
ties arise in finite time in initially smooth solutions to ljl.lt - which motivates one 
to seek for weak solutions understood in the sense of distributions, containing 
for instance shock waves. It turns out that weak solutions are not uniquely 
determined by their initial data at t = 0, say and consequently it is necessary to 
impose an "entropy condition" in order to formulate a well-posed initial value 
problem. 

In the present work, the weak solutions of interest are realizable as limits 
(e — » 0) of smooth solutions if = u e (t, x) to an augmented model of the form 

u\ + f(u% = R$ := (R(e u%, e 2 u\ x , . . .)) x . (1.2) 

The parameter e > introduces a small-scale in the problem, and the term R e x 
accounts for high-order physical features, neglected at the hyperbolic level of 
modeling lll.lb . In the examples arising in fluid dynamics and material science, 
suitable restrictions arise on the right-hand side R e x so that, for instance, the 
conservation laws ( jl.lt are recovered in the limit e — > 0. 

Fine properties of the physical medium under consideration are taken into 
account in l|1.2t . For instance, lll.lt might be the Euler equations for liquid- vapor 
mixtures, while its augmented version dl.2b is a Korteweg-type model includ- 
ing both viscosity and capillarity effects. In many examples of interest, the 
augmented system admits global-in-time smooth solutions that are uniquely 
determined by their initial data. In this context, it is natural to select those 
solutions of (jl.lt that are realizable as (possibly only formal) limits of solutions 
to (|1.2t , that is, to define 

u := \imu £ . (1.3) 

e— »0 
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We face here the fundamental issue raised in the present review: How may 
one characterize this limit u ? In particular, which propagating discontinuities 
should be considered admissible ? Our primary focus is thus on the deriva- 
tion of the "sharp interface" theory associated with l|1.2b consisting of finding 
algebraic conditions for shock waves which, within the class of weak solu- 
tions with bounded variation, single out a unique solution to the initial value 
problem associated with ( jl.H (or, at least, a unique solution to the Riemann 
problem). 

A large part of the mathematical literature on shock waves is focused on 
small amplitude solutions generated by the vanishing viscosity method, corre- 
sponding in l ll.2|l to the regularization 

JS:=ei&. (1.4) 

With this regularization, the limits u := lirrv_>oM e are characterized uniquely 
[154 , 91 , 29 . 8 . 25 1 by one of the equivalent formulations of the "classical" entropy 
condition introduced by Lax (951, Oleinik fl2"9l , Kruzkov [91], Wendroff 11551 , 
Dafermos 15211551 , and Liu |TT9l ll20|. Following |102], we refer to these weak 
solutions as classical entropy solutions. 

1.2 Diffusive-dispersive approximations 

Our objective is to develop a theory of weak solutions that encompasses reg- 
ularizations beyond CE3K including diffusive-dispersive regularizations of the 
form 

K = eu lx + Oie 2 u e xxy .. (1.5) 

Here, a is a fixed real parameter and e — > 0, so that the two regularization terms 
are kept in balance and a subtle competition takes place between two distinct 
phenomena: the diffusion term e u £ xx has a regularizing effect on shock waves 
while the dispersion term a e 2 u e xxx generates oscillations with high frequencies. 
However, as is for instance observed in numerical experiments, such oscillations 
arise near the jump discontinuities, only, and form spikes with finite amplitude 
superimposed on shocks. In turn, the formal limit il.3} does exist but do not 
coincide with the one selected by the viscosity approximation (|1.41 . In fact, 
these limiting solutions depend on the value of the parameter a, and the same 
initial data with different values of a give rise to different shock wave solutions. 

The above observation motivates us to develop, for the system of conser- 
vation laws (Q) , several theories of shock waves, each being associated with 
a different formulation of the entropy condition. In the present work, we are 
primarily interested in hyperbolic systems for which the characteristic speed 
Aj = Aj(u) along each integral curve of the vector field r; = r;(u) admits one 
externum point, at most. For such systems and the class regularization (|1.5t , a 
theory of entropy solution is now available which is based on the concept of a 
kinetic relation. The basic strategy is to impose a single entropy inequality in 
the sense of distributions (cf . | |2.4|| , below), together with an additional algebraic 
condition imposed on certain nonclassical undercompressive shock waves. 



4 



The kinetic relation was first introduced and developed in the context of the 
dynamics of materials undergoing phase transitions, described by the system 
of two conservation laws 

w t -v x = 0, 

i ^ 2 (1-6) 

v t - o(w) x = e v xx - a e w xxx , 

where v denotes the velocity w > — 1 the deformation gradient, and a = a(w) 
the stress of the material. The parameters e and a e 2 represent the (rescaled) 
viscosity and capillarity of the material. The mathematical research on jl.6\ 
began with works by Slemrod in 1984 on self-similar solutions to the Riemann 
problem [ 143 . 144 J (i.e. the initial value problem with piecewise constant data), 
and Shearer in 1986 on the explicit construction of a Riemann solver when a = 
[ 140 [ . The notion of a kinetic relation for subsonic phase boundaries was intro- 
duced by Truskinovsky in 1987 11491 11501 11511 . and Abeyaratne and Knowles 
in 1990 fUGO, wno solved the Riemann problem for il.6i and investigated the 
existence of traveling wave solutions when a is a piecewise linear function. In 
1993, LeFloch [100J introduced a mathematical formulation of the kinetic rela- 
tion for 111. 61 within the setting of functions of bounded variation and tackled 
the initial value problem via the Glimm scheme; therein, the kinetic relation 
was interpreted as an entropy dissipation measure (cf. Section |2j below). The 
kinetic relation was then extended to general hyperbolic systems by the author 
together with many collaborators, and developed into a general mathematical 
theory of nonclassical shocks, which we review in the present notes. For earlier 
reviews, see IToTllT02l . 



1.3 Main objectives of the theory 

Weak solutions in the sense of distributions can be defined within, for instance, 
the space L°° of bounded and measurable functions. However, to deal with 
nonclassical solutions one needs to prescribe (pointwise) algebraic conditions at 
jumps and, consequently, one must work within a space of functions admitting 
traces such as the space of functions with bounded variation. In turn, the 
theory of nonclassical shocks relies on a pointwise formulation of the notion 
of entropy solution and on techniques of pointwise convergence. The main 
difficulty overcome here in our analysis of nonclassical entropy solutions comes 
from their Jack of monotonicity with respect to initial data. The following issues 
have been addressed. 

1. Derivation and analysis of physical models. 

Such nonclassical undercompressive waves arise in weak solutions to a 
large variety of physical models from fluid or solid dynamics: phase tran- 
sitions, thin liquid films, magnetohydrodynamics, Camassa-Holm model, 
martensite-austenite materials, semi-conductors, combustion theory, etc. 

2. Mathematical theory of nonclassical entropy solutions. 
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- By considering initial data consisting of two constant states sepa- 
rated by a single discontinuity, one constructs a nonclassical Rie- 
mann solver associated with a prescribed kinetic function, compati- 
ble with an entropy inequality. 

- Kinetic functions can be determined from traveling wave solutions to 
an augmented model like il.2t - dl.5t and, for certain models, mono- 
tonicity and asymptotic properties of the kinetic functions can be 
established. 

- Establishing via Glimm-type schemes the existence of nonclassical 
entropy solutions to the initial value problem associated with jl-lt 
requires a uniform estimate on the total variation of solutions. Such 
a bound is derived with the help of generalized variation functionals 
that are designed to be diminishing in time. 

- The zero diffusion-dispersion limits for the initial value problem is 
justified rigorously by using weak convergence techniques: compen- 
sated compactness, measure-valued solutions, kinetic formulation. 

3. Approximation of nonclassical entropy solutions. 

- Finite difference schemes with controled dissipation are constructed 
from entropy conservative flux-functions and a careful analysis of 
their equivalent equations. 

- Kinetic functions are associated with finite difference schemes and 
computed numerically. Limiting solutions depend on the regular- 
ization terms arising in the augmented model and, for a given initial 
value problem, different solutions are obtained if different viscosity- 
capillarity ratio or different discretization schemes are used. 

4. Applications and limitations of the theory. 

- Nonclassical solutions exhibit particularly complex wave structures, 
and the kinetic function appears to be the proper tool to represent 
the entire dynamics of nonclassical shocks. 

- Of course, describing all singular limits of a given augmented model 
via a "purely" hyperbolic theory need not always be possible. For 
certain regularizations for which the traveling wave analysis does 
not lead to a unique Riemann solver, a nucleation criterion may 
be required in order to uniquely characterize nonclassical entropy 
solutions. 

One can not underestimate the fruitful connection between the theory of ki- 
netic relations for undercompressive shocks and the so-called DLM (Dal Maso, 
LeFloch, Murat) theory of nonconservative hyperbolic systems. This theory 
was introduced and developed in |98 99. 54] 1 10411 . Fundamental numerical 
issues were first discussed in (761; for recent progress, see (THJ [TH US [35 J and 
the references therein. 
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Building on the pioneering papers ITl ll00[[T49l , the research on undercom- 
pressive shocks developed intensively in the last fifteen years, and the author 
of this review is very grateful to his collaborators, postdocs, and students who 
accompanied him on this subject, including B.T. Hayes (basic theoretical and nu- 
merical issues), N. Bedjaoui (traveling wave solutions), K.T. Joseph (self-similar 
approximations), M.-D. Thanh (Riemann problem), M. Shearer (nucleation cri- 
terion), D. Amadori, P. Baiti, M. Laforest, B. Piccoli (Glimm-type methods), J.C. 
Correia, C. Kondo (vanishing diffusive-dispersive limits), F. Boutin, C. Chalons, 
F. Lagoutiere, J.-M. Mercier, S. Mishra, M. Mohamadian, and C. Rohde (numer- 
ical methods). 

2 Kinetic relations for undercompressive shocks 
2.1 A single entropy inequality 

We follow the strategy advocated in LeFloch [100. 102] for the analysis of the 
formal limits il.3) associated with a given augmented model. Recall that we 
are interested in deriving an entropy condition that singles out all solutions to 
il.li realizable as limits of smooth solutions to {12) . 

In the applications, 111.21 is "compatible" (in a sense defined below) with 
one particular mathematical entropy of the hyperbolic system {1.1) . Hence, we 
assume that {1.1) is endowed with a strictly convex entropy pair, denoted by 
(U,F). (Strict convexity is imposed throughout this presentation, but can be 
relaxed on certain examples.) By definition, (U, F) : 11 — > K X ]R N is a smooth 
map such that 

djF = d k Udjf k , j=l,...,N, 

where dj denotes a partial derivative with respect to the conservative variable 
Uj, with u = (u\, . . .,Un) and / = (f 1 , . . - ,f N ), and we use implicit summation 
over repeated indices. 

It is easily checked from the definition of an entropy that every sufficiently 
regular solution to il.l) satisfies the additional conservation law 

U{u) t + F(u) x = 0. (2.1) 

Definition 2.1. An augmented version (L2) of a system of conservation laws {l.D 
endowed with a strictly convex entropy pair (U,F), is said to be conservative and 
dissipative with respect to the entropy LI if for every non-negative, compactly 
supported, smooth function 6 = 8(t,x) 

lim ff R e 6dtdx = 0, (2.2) 

JJr+xR 

and 

lim ff VU(u £ )-R e 6dtdx<0, (2.3) 
JJr+xr 

respectively. 
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The conditions in Definition 12.11 are fulfilled by many examples of interest 
and, when the singular limit i ll. 31 exists in a strong topology, they lead to the 
conservation laws jl.li and an entropy inequality. 

Definition 2.2. Suppose that il.lty is a system of conservation laws endowed with a 
strictly convex entropy pair (U,F). A bounded and measurable function with locally 
bounded variation u = u(t,x) e U is called an entropy solution if the equations 
together with the entropy inequality 

U(u) t + F(u) x < 0, (2.4) 

hold in the sense of distributions. 

This should be regarded as a preliminary definition of solution since, in 
general, (|2.4l l is not sufficient to guarantee uniqueness for the initial value 
problem, and we are going to impose additional constraints. 

The case of scalar, convex flux-functions is comparatively simpler: one 
entropy inequality turns out to be sufficient to characterize a unique limit 
(see Panov 111321 and De Lellis, Otto, and Westdickenberg 1551 ) and, conse- 
quently, shock waves satisfying a single entropy inequality are regularization- 
independent. The corresponding solutions are referred to as classical entropy 
solutions and contain compressive shocks, only, which satisfy Lax's shock in- 
equalities. 

For strictly hyperbolic systems admitting characteristic fields that are either 
genuinely nonlinear (VAy • r; + 0) or linearly degenerate (VAy • r; = 0), the in- 
equality l|2.4) l for a given strictly convex entropy is sufficiently discriminating to 
select a unique entropy solution. This fact was established first for the Riemann 
problem in Lax's pioneering paper [95]. For the Cauchy problem associated 
with systems with genuinely nonlinear or linearly degenerate characteristic 
fields, a general uniqueness theorem in a class of functions with bounded vari- 
ation was established by Bressan and LeFloch Il29l . Hence, weak solutions 
to genuinely nonlinear or linearly degenerate systems are independent of the 
precise regularization mechanism R £ in the right-hand side of 1 11.21 , as long as 
it is conservative and dissipative in the sense of Definition [2J] 

2.2 Kinetic relation. Formulation based on state values 

In the applications, many models arising in continuum physics (see Section [3J 
below) do not have globally genuinely nonlinear characteristic fields. In numer- 
ical experiments with such systyems, weak solutions often exhibit particularly 
complex wave patterns, including undercompressive shocks. Distinct solutions 
are obtained for the same initial value problem if one changes the diffusion- 
dispersion ratio, the regularization, or the approximation scheme. From the 
analysis standpoint, it turns out that for such systems one entropy inequality 
J2.4I I is not sufficiently discriminating: the initial value problem admits a large 
class of entropy solutions (satisfying a single entropy inequality). Weak solu- 
tions are highly sensitive to small-scales neglected at the hyperbolic level of 
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physical modeling and one needs to determine further admissibility conditions 
beyond | |2.4|| . No universal admissibility criterion is available in this context 
but, instead, several hyperbolic theories must be developed, each being deter- 
mined by specifying a physical regularization. The approach advocated by the 
author is based on imposing a kinetic relation in order to uniquely characterize 
the dynamics of nonclassical undercompressive shocks. 

Recall that a shock wave (u-, u+) is a step function connecting two constant 
states U-, u+ through a single jump discontinuity, propagating at some finite 
speed A = A(m_, u+). To be an entropy solution in the sense of Definition 12. 2[ 
(u_, u+) must satisfy the Rankine-Hugoniot relation 

- A (u + -U-)+ f(u + ) - /(«_) = 0, (2.5) 
as well as the entropy inequality 

-~A~(u(u+) - U(u-j) + F(u + ) - F(u-) < 0. (2.6) 

Denote by £ u the set of all pairs (m_, u+) ellxll satisfying d2.5t -l l2.6l l for some 
speed A = A(u_,w+). 

The emphasis in this section is on solutions to strictly hyperbolic systems 
with small amplitude (but the presentation extends to large data for particular 
examples). Then, one can show that, for each U-, 02.51 consists of N curves 
issuing from U- and tangent to each characteristic vector ry(w_). Wen the shock 

speed A = A(u-, u+) is comparable to the characteristic speeds Aj(u±), we write 
also A = Aj(u-, u+) and accordingly we decompose the shock set 



Furthermore, we tacitly assume that the set of definition IX is replaced by a 
smaller open subset, if necessary. Finally, we assume that the characteristic 
speed Aj = A;(u) along each integral curve of the vector field r; = rdu), or 
along each Hugoniot curve, admits one extremum point, at most. We refer to 
characteristic fields having such an extremum as concave-convex characteristic 
fields. 

A /-shock (w_, u+) is called slow undercompressive if 

Aj(u±) > Aj(u-,u+) 

and fast undercompressive if the opposite inequalities hold. Across under- 
compressive shocks, one must supplement the Rankine-Hugoniot relation and 
entropy inequality with an additional jump condition, as follows. Recall that 
functions of bounded variation (BV) admit traces in a measure-theoretic sense. 
While 02.11 and 02.41 can be imposed in the sense of distributions, this regularity 
of BV functions is here required for pointwise conditions to make sense. 

Definition 2.3 (Formulation based on state values). Consider small amplitude so- 
lutions to a strictly hyperbolic system of conservation laws {1.1) endowed with a strictly 
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convex entropy pair (U,F) and admitting genuinely nonlinear, linearly degenerate, or 
concave-convex characteristic fields, only. 

1. A kinetic function compatible with the entropy U for a concave-convex 
j-characteristic family is a Lipschitz continuous map <p. = qr.{u) defined over M such 
that 

(u,<pXu)) e is an undercompressive shock. (2.7) 

2. A family of kinetic functions qft := (<p^) being prescribed for each concave- 
convex characteristic field, one says that a bounded and measurable function with 
locally bounded variation u = u(t,x) eliisa tp^-admissible entropy solution if it is 
an entropy solution to l|l.lb - (|2.4l l and, at every point of approximate jump discontinuity 
(t,x) ofu associated with an undercompressive j-shock, 

u + = cp){u-), (2.8) 

where U-, u + are the left- and right-hand limits at that point. 

We will refer to (|2.8b as the kinetic relation associated with cp^, and to in- 
admissible entropy solutions as nonclassical entropy solutions. Note that not 
all propagating waves within a solution require a kinetic relation, but only 
undercompressive shocks do. The kinetic relation is very effective and, for 
many (but not all) models of interest, selects a unique cp^-admissible entropy 
solution to the Riemann problem. 

2.3 Kinetic function. Formulation based on the entropy dissi- 
pation 

To have an effective theory, the kinetic function must be specified. One may 
simply postulate the existence of the kinetic function (p° without refereeing to 
a small-scale modeling, but instead defining it via laboratory experiments or 
some physical heuristics. A sounder approach relies on a prescribed augmented 
model, as we now explain. 

We introduce the kinetic function in the following way based on the entropy 
dissipation measure generated by a given augmented model. Suppose that the 
product VU(u e ) ■ R e arising in the right-hand side of \1.2\ can be decomposed 
in the form 

VU(u e ) -R £ = Q £ - p e , 

where Q £ converges to zero in the sense of distributions and p e is a uniformly 
bounded sequence of non-negative L 1 functions. Let us refer to \f as the entropy 
dissipation measure. Note that it depends on the specific regularization R e and 
the entropy U. (This decomposition can be established for several examples; 
see Section [3]) After extracting a subsequence if necessary, this sequence of 
measures converges in the weak-star sense 

p u (u) := limp e < 0, 
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where the limiting measure p. u (u) is clearly related to the pointwise limit 
u := lim e _>oM £ / but in general cannot be uniquely determined from the sole 
knowledge of this limit. 

For regularization-independent shock waves the sole sign of the entropy 
dissipation measure fi u (u) suffices and one simply writes down the entropy 
inequality J2.4I I. However, for regularization-sensitive shock waves, the range 
of the measure p u (u) plays a crucial role in selecting weak solutions. It is 
proposed in LeFloch H100U101l [T02| to replace the entropy inequality d2.4b by 
the entropy equality 

U(u) t + F(u) x = jj u (u) < 0, (2.9) 

where fi u (u) is a non-positive, locally bounded measure depending on the 
solution u under consideration. Clearly, the measure p u (u) cannot be prescribed 
arbitrarily and, in particular, must vanish on the set of continuity points of u. 

Recalling that we focus the presentation on strictly hyperbolic systems and 
small data, let Aj be the range of the speed function A/, and let us use the 

short-hand notation £ u X A for the union of all subsets £. u X A,. 

1 ' 

Definition 2.4 (Formulation based on the entropy dissipation measure). Sup- 
pose that (1J) z's a system of conservation laws endowed with an entropy pair (U,F). 

1. A kinetic function (compatible with the entropy LI) is a Lipschitz continuous 
map O = 0(w_, u+, A) defined on £ u x A and satisfying 

O < 0. (2.10) 

2. Let O be a kinetic function. A bounded and measurable function with locally 
bounded variation u = u(t, x) e U is called a <5-admissible entropy solution if it is an 
entropy solution of jl.l^ - HAi and if at every point of approximate jump discontinuity 
(t,x) ofu 

- A (U(u+) - U{u-)) + F(u+) - F(u-) = 0(m_, u+, A), (2.11) 

where m_, m+, A denote the left- and right-hand limits and A the shock speed. 

The entropy dissipation (|2.11l l is analogous to what is called a "driving force" 
in the physical literature. Note that at a point (w_, u+) where the following trivial 
choice is made 

0(m_, u+, A) := -A (U(u + ) - ii(u_)) + F(u + ) - F(w_), (2.12) 

the condition Il2.11b is vacuous. In practice, 0(m-, u+, A) will be strictly monotone 
in A for all undercompressive shocks, and Il2.11b may be seen to select a unique 
propagation speed A for each m_. (See [102[ for details.) 

To determine the kinetic function, one then analyzes traveling wave solu- 
tions associated with \1.2\ . Given any shock wave («_, u+), a function u e (t, x) = 
w{y) with y := (x — A t)/e is called a traveling wave associated with the shock 
(u_, u+) if it is a smooth solution to (|1.2) l satisfying 

w(-oo) = «_, w(+oo) = u + , (2.13) 
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and 

lim w'{y) = lim w"(y) = . . . = 0. (2.14) 

|y|— >+oo |y|— »+co 

The function w satisfies the following system of ordinary differential equations 

R{w', w", ...) = -A(w - u-) + f(w) - f(u-). (2.15) 

At this juncture, it is straightforward, but fundamental, to observe the following 
properties. 

Proposition 2.5 (Kinetic relations derived from traveling waves). Consider an 
augmented version jl.2i of a system of conservation laws 1 11. 11 which is endowed with 
a strictly convex entropy pair (U, F), and is conservative and dissipative with respect to 
U. Then, given a traveling wave solution satisfying <|2.13b — <|2.15b , the pointwise limit 



u(t,x) := \imw(- — — ) = { 



e^O V £ ' \U+, X > A t, 

is a weak solution to jl.li satisfying the entropy inequality il.ty . In particular, the 
Rankine-Hugoniot relation H2.5I I follows by letting y — > +oo. 

Moreover, the solution u satisfies the kinetic relation j2.11|l in which the dissipation 
measure reads 

p u (u)=M u [u]5 x - At , 

M u [u] ■=- \ R(w(y), w'{y), w"(y), . . .) • V 2 U(w) w'(y) dy, 
Jk 

where b x -\ t denotes the Dirac measure concentrated on the line x — A t — 0. 

3 Physical models 

3.1 Nonlinear diffusion model 

We begin with the case of a conservation law regularized with diffusion, i.e. 

u\ + f{u% = e (b(u e ) ul) , u e = u e (t, x) e U, (3.1) 

where the function b : R — > (0, oo) is bounded above and below, and e > is 
a small parameter. The following observation goes back to Kruzkov (9T1l and 
Volpert II 15411 and, in particular, given suitable initial data, the solutions to the 
initial value problem associated with [ 13.1| | converge strongly as e — > 0. 

Lemma 3.1 (Nonlinear diffusion model). Solutions u e to the augmented model 
J3.1b satisfy, for every convex function U : ]R — > ]R, 

U(u £ ) t +F(u e ) x = -D £ + eC x . 

D e := e b(u £ ) U"(u e ) \u £ x \ 2 , C := b(u £ ) U(u £ ) x , 
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in which F(u) := J" f'(v)U'(v)dv. Hence, the limit u = \im e ^ou e satisfies the 
following entropy inequalities associated with the nonlinear diffusion model 

U(u) t + F(u) x < 0. (3.2) 

We refer to weak solutions u e L°° to the hyperbolic conservation law 

u t + f(u) x = 0, u = u(t, x)eR (3.3) 

that satisfy all of the inequalities i3.2) as classical entropy solutions. Clearly, 
the inequalities (|3.2|l are equivalent to the Kruzkov inequalities 

\u - k\ t + (sgn( H - k)(f(u) - f(k))) x < 0, k 6 R 

Provided additional regularity beyond L°° is assumed and, for instance, when 
u has bounded variation, the classical entropy solutions to ( 13.31 1 can be also 
characterized by the Oleinik inequalities (cf. 14.21 below), which, at shocks, 
impose a local concavity or convexity property of /. 

3.2 Linear diffusion-dispersion model 

Nonclassical solutions are obtained when both diffusive and dispersive effects 
are taken into account. A typical model of interest is provided by the following 
conservation law 

«t + f(u% = e 4 X + 7(e) u £ xxx , u e = u e (t, x), (3.4) 

in which e > and y — y{e) are real parameters tending to zero. This equation 
was studied first by Shearer et al. |85L Hayes and LeFloch 1 71 1, and Bedjaoui 
and LeFloch EH . 

The relative scaling between e and y(e) determines the limiting behavior of 
the regularized solutions, as follows: 

• When y(e) « e 2 , the diffusion plays a dominant role and the limit u := 
lim e _>o u e is a classical entropy solution. 

• When y(e) » e 2 , the dispersion effects are dominant and high oscillations 
develop in the limit, allowing only for weak convergence of u £ . The 
theory for dispersive equations developed by Lax and Levermore |97] is 
the relevant theory in this regime. 

• Finally, in the balanced regime where y(e) := oce 2 for fixed a, the limit 
u :- lirrv^o u e exists in a strong sense and solely mild oscillations arise 
near shocks. The limit u is a weak solution to the hyperbolic conservation 
law d3.3l >. When a > 0, the limit u exhibits a nonclassical behavior, and 
strongly depends on a. 
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The above observations motivate us in the rest of this paper, and we adopt 
the following scaling 

u e t + f(u e ) x = eu e xx + ae 2 u e xxx , u e = if = u e (t,x), (3.5) 

where a is fixed. Following Hayes and LeFloch [71 |, we now derive one entropy 
inequality in the limit. The contribution due to the diffusion decomposes into a 
non-positive term and a conservative one; the contribution due to the dispersion 
is entirely conservative; hence, formally at least, as e — > we recover the entropy 
inequality. 

Lemma 3.2 (Linear diffusion-dispersion model). For the augmented model \3A\ , 
one has 

(l/2)\u 2 e \ t + F(u% = -D £ + eC x , 

D e := e \u e x \ 2 > 0, C := u e u% + a e (u e u e xx - (1/2) \u x \ 2 ). 

Hence, as e — > 0, the (possibly formal) limit u := lim^o u e satisfies the single entropy 
inequality associated with the linear diffusion-dispersion model 

(u 2 /2j t + F(u) x < 0, F':=uf. 

No specific sign is available for arbitrary convex entropies. 

More generally, consider the nonlinear diffusion-dispersion model 

K? + f{u e ) x = e(b(u £ ) u e x ) x + ae 2 (c a ( M e ) (c 2 {if) u e x )^, (3.6) 

where the functions b, C\,Ci are given smooth and positive functions, and a is a 
real parameter. 

Lemma 3.3 (Nonlinear diffusion-dispersion model). For the augmented model 
l|3.6b , the formal limit u = lim e _,o u e satisfies the following entropy inequality asso- 
ciated with the nonlinear diffusion-dispersion model 

U(u) t + F{u) x < 0, 

U" = £i>0, F':=f'W. 

c\ 

Proof. In the entropy variable u = U'(u), the dispersive term takes the symmet- 
ric form 

(ci(«) (c 2 (w) u x j J x = \ci{u){ci{u)Hx)J x 
and, consequently, any solution of d3.6b satisfies 

U{u% + F(u e ) x = -D e + eC x . 
D e :=eb(u) U"(u)\u x \ 2 , 

C := b(u) U'{u)u x + ae(ci{u)u(ci{u) u x ) - \c 2 {u) u x \ 2 /2j. 

□ 
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3.3 Thin liquid film model 

Consider next the augmented model 

u\ + (u 2 e - u% = e (ulu e x ) x - y(e) (u 3 u e xxx ) x (3.7) 

with e,y = y(e) > 0, in which the right-hand side describes the effects of 
surface tension on a thin liquid film moving on a surface. In this context, 
u = u(t,x) e [0, 1] denotes the normalized thickness of the thin film layer, while 
parameters governing the various forces and the slope of the surface are typi- 
cally incorporated into the parameter e. More precisely, the equation ( 13.71 1 arises 
from the lubrication aproximation of the Navier-Stokes equation and models 
the physical situation in which the film is driven by two counteracting forces: on 
one hand, the gravity is responsible for pulling the film down an inclined plane 
and a thermal gradient (resulting from the surface tension gradient) pushing 
the film up the plane. The thin liquid film model was studied by Bertozzi and 
Shearer (22), together with Munch (201, Levy 0115111161 , and Zumbrun |2l] See 
also Otto and Westdickenberg 111301 , and LeFloch and Mohamadian [107]. 

For the purpose of the present paper we need the following observation, 
made in ITTTUI . 

Lemma 3.4 (Thin liquid film model). For the augmented model 113.71 1, one has 

(if log u £ - u e ) t + ((u 2 e - ul) log u e - u e + u 2 .) = -D e + e C x 

with D e := e u\ |m^| 2 + y(e) |(m 2 u x ) x \ 2 > 0. In the limit e — > one deduces the single 
entropy inequality associated with the thin liquid film model 

(u log u - u) t + ((w 2 - m 3 ) log u-u + uj < 0. 

3.4 Generalized Camassa-Holm model 

Consider now the following conservation law 

u e t + f(u e ) x = e u xx + a e 2 (u e txx + 2u £ x u xx + u e u xxx ), (3.8) 

in which a is a fixed parameter and e — > 0. This equation arises as a simplified 
model for shallow water when wave breaking takes place, and was studied 
by Bressan and Constantin [28 1, Coclite and Karlsen |40], and LeFloch and 
Mohamadian 11071 . 

Lemma 3.5 (Generalized Camassa-Holm model). For the augmented model <|3.8b 
one has 

((\if\ 2 + ae 2 K| 2 )/2), + F(u% = -e \uf + e C x , 

which in the limit e — » implies the following entropy inequality associated with 
the generalized Camassa-Holm model 

(w 2 /2) t +F( M ). Y <0, F':=uf. 



15 



The entropy inequality coincides with the one in Lemma l3T2l for the diffusion- 
dispersion model. As pointed out in [107[, limiting solutions look quite similar 
but still do not coincide with the ones obtained with the diffusion-dispersion 
model. 



3.5 Phase transition model 

Let us return to the model 1 11. 6b cited in the introduction which, without regu- 
larization, reads 

w t - v x = 0. 

< \ n (3 ' 9) 
Vt - a{w) x = 0. 

For many typical elastic materials, we have 

a' (to) > for all w > -1 (3.10) 

so that l|3.9|l is strictly hyperbolic with two distinct wave speeds, — Aj = A2 = 
c(w) (the sound speed). The two characteristic fields are genuinely nonlinear 
if and only if a" never vanishes. However, many materials encountered in 
applications do not satisfy this condition, but rather loose convexity at w = 0, 
that is, 

a"(w) % if w ^ 0. 

One mathematical entropy pair of particular interest is the one associated with 
the total energy of the system and given by 

v 2 

U(v, w) =— + L(w), F(v, w) = —o(w) v, 

1 (3.11) 



-J 

Jo 



H(w) := I o(s)ds. 



The entropy U is strictly convex under the assumption (|3.10l l. 

Material undergoing phase transitions may be described by the model \3.9\ 
but with a non-monotone stress-strain function satisfying 

o'(w) > 0, we (-1, w m ) U {w M , +00), 

(3 12) 

o'(w)<Q, zve(zv m ,zv M ) 

for some constants w™ < w M . In the so-called unstable phase (w m ,w M ) the 
system admits two complex conjugate eigenvalues and is elliptic in nature. 
However, the solutions of interest from the standpoint of the hyperbolic theory 
lie outside the unstable region. The system is hyperbolic in the non-connected 
set U:=(Rx (-l,io m )) u(Kx (h^+oo)). One important difference with the 
hyperbolic regime is about the total mechanical energy 03.111 1, which is still 
convex in each hyperbolic region, but (any extension) is not globally convex in 
(the convex closure of) U. 
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The system (|3.9H3.12| | an d its augmented version (|1.6b leads to complex 
wave dynamics, including hysteresis behavior, and is relevant to describe phase 
transitions in many different applications involving solid-solid interfaces or 
fluid-gas mixtures. 

Lemma 3.6 (Phase transition model). For the augmented model jl-61 , one has 

= e (v v x j - e v\ + a e 2 (v x w x - v w xx ^ , 

so that in the limit one formally obtains the following entropy inequality associated 
with the phase transition model 

2 

[-+Z(w)) r (va(w)) x <0. (3.13) 
3.6 Nonlinear phase transition model 

More generally, assume now an internal energy function e = e(w, w x ), and let 
us derive the field equations from the action 



7(y) 



I I (e(zv, w x ) - — jdxdt. 
Jo Jq 2 



Precisely, considering the unknown functions V and W and defining the total 
stress as 

_ de i de \ 

£(w, u? x , w„) : = —( w ,w x ) - [—(w,w x )J^, 

we can obtain 

Vt - L(w, W x , Wxx)x - 0, 
Wt - Vx - 0. 

Including next a nonlinear viscosity /.( = f.i(zu), we arrive at the nonlinear phase 
transition model which includes viscosity and capillarity effects: 

Wt - v x = 0, 

vt - Z{w, w x , zo X x)x = [n(ro)v x ) . 

Again, the total energy E(w,v,w x ) := e{w,w x ) + v 2 /2 plays the role of a 
mathematical entropy, and we find 

E(W, V, W x )t ~ W x , W xx ) p)^ = (v x ^-( w > w >)) x + (f i w ) v v x) x ~ P( w ) v l> 

and once more, a single entropy inequality is obtained. 
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We now specialize this discussion with the important case that e is quadratic 
in w x . (Linear term should not appear because of the natural invariance of the 
energy via the transformation x i-» -x.) Setting, for some positive capillarity 
coefficient A(w), 



2 



e(w,w x ) = e(w) + A(w) 



the total stress decomposes as follows: 

L(w, w x , w xx ) = a(w) + A'(w) — (A(w) w x ) x , a(w) = e'(w), 

and the field equations take the form 

Wt - v x = 0, 

v t - o{w) x = (A (w) — - (A(w) Wx)J x + (w) v x ) x . 

Lemma 3.7 (Nonlinear phase transition model). For the augmented model d3.14t , 
one finds 

(e(w) + y + A(w) ^) ( - (a(w) v) % 

= (/.i(ro) v v x ^ - ji(w) v\ + {v - jj^- w\ - v (A{w) Wxj^ + v x A(iv) w x ^j , 

which leads to the same entropy inequality for the nonlinear phase transition 
model as | |3.13l l in Lemma \3£\ 

When the viscosity and capillarity are taken to be constants, we recover 
the example in Section 1331 above. The entropy inequality is identical for both 
regularizations. 



3.7 Magnetohydrodynamic model 

Consider next the following simplified version of the equations of ideal mag- 
netohydrodynamics 

Vt + ((v 2 + w 2 ) V) = £ v xx + a £ zv xx , 

f,2 7 \ (315) 

w t + y(v + w ) u>J=e w xx - a e v xx , 

where v, iv denote the transverse components of the magnetic field, e the mag- 
netic resistivity, and a the so-called Hall parameter. The Hall effect taken into 
account in this model is relevant to investigate, for instance, the Earth's solar 
wind. 

When a = 0, 1 13.151 was studied by Brio and Hunter [30], Freistuhler et al. 
|62j|63j|6!), Panov 11311 , and others. The equations are not strictly hyperbolic, 
and for certain initial data the Riemann problem may admit up to two solutions. 
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When a + 0, we refer to LeFloch and Mishra [106] who demonstrated 
numerically that a kinetic function can be associated to this model. For the 
purpose of the present section, we observe the following. 

Lemma 3.8 (Magnetohydrodynamic model). For the augmented model J3.15I I one 
has 

(1/2) (v 2 £ + wl) t + (3/4) ((v 2 e + wl)\ = - e ({vl) 2 + K) 2 ) + e C%, 

so that in the limit e — > the following entropy inequality associated with the 
magnetohydrodynamic model holds 

(1/2) (v 2 + w 2 ) t + (3/4) {(v 2 + w 2 ) 2 )^ < 0. 



3.8 Other physical models 

The Buckley-Leverett equation for two-phase flows in porous media provides 
another example, studied in Hayes and Shearer [74J and Van Duijn, Peletier, and 
Pop 11531 . Finally, we list here several other models of physical interest which, 
however, have not yet received as much attention as the models presented so far. 
Since the hyperbolic flux part of these models admits an inflection point and, 
moreover, physical modeling includes dispersive-type terms, it is expected 
that compressive and undercompressive shocks occur in weak solutions, at 
least in certain regimes of applications. The actual occurrence of nonclassical 
shocks depends upon the form of the regularization, and further investigations 
are necessary about the quantum hydrodynamics models |H l86l , phase field 
models I3TII134 I, Suliciu-type models |32, 65, 145], non-local models involving 
fractional integrals [89, 135, 136], and discrete molecular models based, for 
instance, on potentials of the Lennard-Jones type l23ll57lll28llT52ll59l . 



4 Nonclassical Riemann solver 



4.1 Consequences of a single entropy inequality 

For simplicity in the presentation, we consider a scalar equation 

Ut + f(n) x = 0, u = u(t,x) (4.1) 
with concave-convex flux 

uf"(u)> forw^o, 
f"'(0) * 0, lim f'(u) = +co. 



It— >±oo 



To the flux /, we associate the tangent function qfi : ]R — » ]R (and its inverse 
denoted by <p~^) defined by 

, h f(u) - f((p\u)) 



19 



Weak solutions u 6 L°° to (|4.11 , by definition, satisfy 
^J" <pf + /(u) (px^j dxdt = 



for every smooth, compactly supported function (p. If u is a function with 
locally bounded variation, then Ut and u x are locally bounded measures and 
||4.1| | holds as an equality between measures. A shock wave (m_, u+) e R 2 , given 
by 

\u-, x < At, 
u(t,x) = i 

\u+, x> At, 

is a weak solution to (|4.1b provided the Rankine-Hugoniot relation 

-A(w + -w_)+/(w + )-/(w_) = 

holds. In the scalar case, this relation determines the shock speed uniquely: 

/(«_)- /( M+ ) 

A = =: a(u-,u+). 

U- - u+ 

Motivated by the physical models studied in the previous section, we im- 
pose that weak solutions satisfy a single entropy inequality 

U(u) t + F(u) x < 0, U" > 0, F» := f'(u) U'(u) 

for one prescribed entropy pair. In other words, on the discontinuity (u_, u+) 
we impose 



f(u-) - f(u+) I \ 
E(u_, u + ) ■- ~ _ (U(u + ) - U(u-)) + F(u + ) - F(m_) 

<0. 



It is easily checked that 

v< ^ P,r", w J f(v)-f(u-) f(u + )-f(u.) \ A 

E(u-,u+) = - !J (p) (o — m_) I hfo. 

J„ \ v - u- u+-u- ) 

Note in passing that imposing all of the entropy inequalities would lead to 
the Oleinik's entropy inequalities 

m-f(u+) < f(u + )-f(u-) 

v - u+ u+ - li- 

fer all v between U- and u+. Imposing a single entropy inequality is much 
weaker than imposing (|4.2t . 
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Proposition 4.1 (Zero entropy dissipation function). Given a concave flux f and 
a strictly convex entropy U, there exists a function qr Q : M i-» R such that 

e(u, <Po(«)J = 0, <p b (u) + u ( when u + 0) 
(<Po ° <Po)( M ) = 

Moreover, for instance when u > 0, f/ie entropy dissipation E\u,q>g(u)j is negative if 
and only ifue {cp^{u), uj. 

The above result follows from the identity (for u_ + u+): 

d u+ E(u-,u+) = b(u-,u+)d u+ a(u-,u+), 

b(u-, u+) := U(u-) - U(u+) - U'(u+) (u- - u+) > 0, 

where the sign of the factor 

-, , f'(u + )-a(u-,u+) 

d u+ a(u-,u+) = 

u+ - it- 
is easily determined in view of the concave-convex shape of /. 

4.2 Admissible waves 

We are in a position to deal with the Riemann problem, corresponding to the 
initial data 

Ui, x < 0, 



u(x, 0) = 



u r , x > 0, 



with (w;, Uy) e R 2 . It turns out that a single entropy inequality allows for three 
types of waves (m_, u+): 

• Classical compressive shocks, having 

u- > 0, (p\u~) < u+ < U-, 

which satisfy Lax shock inequalities 

f(u+) - f(u-) 
f ^>~ u+ - u_ ' 

Compressive shocks arise from smooth initial data: for instance, using 
the method of characteristics one can write 



i{t,x) = u(0,x-tf'(u(t,x))), 



and one sees that the implicit function theorem may fail to determine 
the value u(t, x) uniquely, whenever t is sufficiently large. Compressive 
shocks arise also from singular limits, for instance from vanishing viscos- 
ity limits. 
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Nonclassical undercompressive shocks, having 

u- > 0, <Pq(m-) < u+ < <p\u-), 
for which all characteristics pass through it: 

/(«+)-/(«-) 



min(/'( M _),/'( M+ ))> 



u+ - U- 



The cord connecting u_ to u+ intersects the graph of/. Undercompressive 
shocks arise from certain (dispersive) singular limits, only. Il3.5t . 

• Rarefaction waves, which are Lipschitz continuous solutions u depending 
only upon E, :— x/t and satisfy the ordinary differential equation 

+/(m(£)) c =0. 

Precisely, a rarefaction consists of two constant states separated by a self- 
similar solution: 

(U_, X < tf'(U-), 

ifTHx/t), t /'(«_) < x < t f{u-), 
U+, X > t /'(«+). 

This construction makes sense provided /'(«-) < /'(«+) and /' is strictly 
monotone on the interval limited by w_ and u+. 

By attempting to build a solution to the Riemann problem that uses only 
the above admissible waves, one realizes that the Riemann problem may ad- 
mit a one-parameter family of solutions satisfying a single entropy inequality. 
Indeed, within an open range of initial data, one can combine together an 
arbitrary nonclassical shock plus a classical shock. 

4.3 Entropy-compatible kinetic functions 

At this stage of the discussion we introduce an additional admissibility require- 
ment. 

Definition 4.2. A kinetic function is a monotone decreasing, Lipschitz continuous 
function <p" : R i— > ]R satisfying 

cpliu) < (p\u) < (p\u), u > 0. (4.4) 

The kinetic relation 

w+ = <p b («-) 

then singles out one nonclassical shock for each left-hand state 

Equivalently one could prescribe the entropy dissipation rate across under- 
compressive shocks. The following two extremal choices could be considered: 
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• In the case <p" = qfi, the kinetic relation selects classical entropy solutions 
only, which in fact satisfy all convex entropy inequalities. 

• The choice (p° = </9g is not quite allowed in l|4.4l l and would correspond to 
selecting dissipation-free shocks, satisfying an entropy equality. 

The property (cp b Q o cp b Q )(u) = u (see (|4.3l l) implies the contraction property 

\<p b (<p\u))\ < \u\, u + 

which turns out to prevent oscillations with arbitrary large frequencies in solu- 
tions to the initial value problem. Introduce the companion threshold function 
qr : 1R — > ]R associated with <p° , defined by 

f(u)-f(cp\u)) f(u)-f(cpt(uj) 

I = 1 ~> u *°- 

u - (p b (u) u - <p*{u) 

Definition 4.3. The nonclassical Riemann solver associates to any Riemann data 
Ui, u r the following entropy solution (for u\ > 0, say): 

• a rarefaction wave ifu r > u\, 

• a classical shock ifu r e [(p"(u;), uA 

• if u r e ^(p b (ui),(p^(uifj r a nonclassical shock (ui,(p b (ui) J followed by a classical 
shock (cp^(ui), u r ^j, and 

• if u r < <p b (ui), a nonclassical shock (ui, <p b (ui)j followed by a rarefaction wave 
(<p b (ui),u r ). 

In conclusion, given a kinetic function (p° compatible with an entropy, the 
Riemann problem admits a unique solution satisfying the hyperbolic conserva- 
tion law, the Riemann initial data, the single entropy inequality, and a kinetic 
relation u + = qr{ii-). 

Observe the L 1 continuous dependence property satisfied by any two en- 
tropy solutions u, v associated with the same kinetic function, 

IK*) - v(t)\\ U{K) < C(T,K) IHO) - v(0)\\ L1{K) 

for all t e [0, T] and all compact set KcR. However, no pointwise version of 
this continuous dependence property holds, as the Riemann solution contain 
"spikes": some intermediate states depend discontinuously upon the initial 
data. 
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4.4 Generalization to systems 

Nonclassical Riemann solvers are also known for several systems of interest. 
The 2x2 isentropic Euler, nonlinear elasticity, and phase transition systems were 
studied extensively in the mathematical literature. The nonclassical Riemann 
solver was constructed by Shearer et al. [139. 142] (for the cubic equation) and 
LeFloch and Thanh lllll[|112[|113 il (for general constitutive equations admitting 
one inflection point). Uniqueness is also known for the Riemann problem 
(within the class of piecewise smooth solutions) when this system is strictly 
hyperbolic. However, in the hyperbolic-elliptic regime, two solutions are still 
available after imposing the kinetic relation. About the qualitative properties of 
solutions, see also important contributions by Hattori [68 , 69 , 70], Mercier and 
Piccoli CUES, and Corli and Tou geron |48j. Partial results are also available 
for the 3x3 Euler equations for van der Waals fluids 11141 . Finally, for the 
construction of the nonclassical Riemann solver to general strictly hyperbolic 
systems of N > 1 conservation laws, we refer to Hayes and LeFloch [73[. 

5 Kinetic relations associated with traveling waves 
5.1 Traveling wave problem 

It was explained in the previous section that kinetic functions characterize the 
dynamics of nonclassical shocks and allow one to solve the Riemann problem. 
The actual derivation of a kinetic relation is an essential issue and, in the present 
section, we explain how to derive it effectively from an analysis of traveling 
wave solutions to a given augmented model. In some cases, the kinetic function 
is computable by (explicit or implicit) analytic formulas. 

For simplicity in the presentation, we consider the augmented model 

u t + f{u) x = a (\u x \ p u x j x + Uxxx (5.1) 

in which / is assumed to be a concave-convex function while a > and p > 
are prescribed parameters. 

Searching for solutions u(t, x) = w{y) depending only on the variable y = 
x — At, we arrive at the second-order ordinary differential equation 

- A (w - «_) + f(w) - f(u-) = a \w'f w' + w" (5.2) 

with boundary conditions 

lim w(y) = u±. 

y->±oo 

Here, u ± and A are constant states satisfying the Rankine-Hugoniot relation. 

Remark 5.1. 1. The parameter e has been removed by reseating of the original 
augmented equation. 

2. All of the results in this section extend to the more general model 

u t + f(u) x = a (b(u, u x ) \u x \ p u x J x + (ci(w) (c 2 {u) u x JJ x , 
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where the functions b(u, v), C\{u), €2(11) are continuous and positive, and b(u, v) \vfv is 
monotone increasing in v. 

In view of J5.2|l , the fundamental questions of interest are the following 
ones: do there exist traveling wave solutions associated with classical and/or 
with nonclassical shock waves («_, «+) ? Can one associate a kinetic function 
(p\ „ to this model ? If so, is this kinetic function monotone ? What is the local 

behavior of <p b at u = ? How does (p^ depend upon the two parameters a, p ? 

Answers to these questions were obtained first for the cubic flux function, 
by deriving explicit formulas for the kinetic function: cf . for p = 0, Shearer et al. 
[85J and, for p = 1, Hayes and LeFloch l7ll . In the latter case, it was observed 
that <p^;i(0) = <Pg (0) = -1. General flux-functions and general regularization 
were covered by Bedjaoui and LeFloch in the series of papers [11J-[15|. 

5.2 Existence and asymptotic properties of kinetic functions 

To show the existence of the kinetic function, we reformulate ||5.2| | as a first- 
order system in the plane (w, w'). The left-hand state and the speed being fixed, 
the corresponding equilibria are the solutions U\ < U2 < U3 = U- to 

-A (w - u-) + f{w) - f(u-) = 0, 

which admits two non-trivial solutions (beyond «_). Equilibria may be saddle 
points (two real eigenvalues with opposite signs) or nodes (two eigenvalues 
with same sign). A phase plane analysis shows that there exist saddle-node 
connections from U- to U2 (corresponding to classical shocks) as well as saddle- 
saddle connections from w_ to Mi (corresponding to nonclassical shocks). To 
state the results precisely, we introduce the following set of all admissible shocks 

S(m_) := |m+ / there exists a TW connecting u± J 

The following theorem is established in Bedjaoui and LeFloch [14J, and shows 
that to the augmented model one can associate a unique kinetic function which, 
furthermore, is monotone and satisfies all the assumptions required in the 
theory of the Riemann problem. 

Theorem 5.2 (Existence of the kinetic function). For each a > and p > 0, 
consider the traveling wave problem j5.2i for the augmented model j5.1) . Then, there 
exists a kinetic function cp b ap : ]R — > R which is locally Lipschitz continuous, strictly 
decreasing, and such that, for instance when u > 0, 

S(u) = {(pl p (u)}u((pi p (u),u], 

<Po(") < <Pip( u ) ^ (p\u), 
where is the companion function associated with cp^ p . Moreover, for < p < 1/3, 
there exists a threshold function Ap satisfying 

Ap : ]R — > [0, 00) Lipschitz continuous, A p (0) = 0, 
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so that 

<p^(u) = qfi{u) if and only if a > Aj,(u), 

implying that all shocks of sufficiently small strength are classical. On the other hand, 
for for all p > 1/3, one has 

<p\u) * (p\u) (u * 0), 
implying that there exist nonclassical shocks of arbitrarily small strength. 

The following local behavior is relevant in the existence theory for the initial 
value problem. 

Theorem 5.3 (Asymptotic properties of the kinetic function). Under the assump- 
tions and notation of Theorem 15.21 the behavior of infinitesimally small shocks is 
described as follows: 

• p = 0: 

<p b '(0) = ^'(0) = -1/2, 
Aj(0) = 0, Aj'(0±) * 0; 

• < p < 1/3: 

/(0) = -1/2, 

A>) = 0, A}'(0±) = +«>; 

• l/3<p<l/2: 

</(0) = -1/2; 

• p = 1/2: 

<p b '(0) G (<p- b '(0), -1/2) = (-1,-1/2), 
lim <p b '(0) = -1, lim <p b '(0) = -1/2; 

a— >0+ a— >+«j 

• p > 1/2: 

/(0) = -1. 

Remark 5.4. Explicit formulas for the kinetic function 185. 71\\14h are available when 
the flux is a cubic, say f(u) = u 3 and p = 0, 1/2 or 1. In particular, when p = 1/2 the 
kinetic function turns out to be the linear function 

<p\u) = -c a u, c a e (1/2,1). 
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5.3 Generalization to systems 



The existence and properties of traveling waves for the nonlinear elasticity 
and the Euler equations are known in both the hyperbolic 111391 H2l and the 
hyperbolic-elliptic regimes 1 150 . 13l [16lH42l . For all other models, only partial 
results on traveling waves are available. 

The existence of nonclassical traveling wave solutions for the thin liquid 
film model is proven by Bertozzi and Shearer in 11221 . For this model, no 
qualitative information on the properties of these traveling waves is known, 
and, in particular, the existence of the kinetic relation has not been rigorously 
established yet. The kinetic function was computed numerically in LeFloch 
and Mohamadian [107J. For the 3x3 Euler equations, we refer to 1 151 . 

The Van de Waals model admits two inflection points and leads to multiple 
traveling wave solutions. Although the physical significance of the "second" 
inflection point is questionable, given that this model is extensively used in 
the applications it is important to investigate whether additional features arise. 
Indeed, it was proven by Bedjaoui, Chalons, Coquel, and LeFloch 1 10J that non- 
monotone nonclassical traveling wave profiles exist, and that a single kinetic 
function is not sufficient to single out the physically relevant solutions. 

The Van der Waals-type model with viscosity and capillarity included reads 

T f - U x = 0, 

«t + V( T )x = «(jS(t) |t x | ? u x } x - t xxx , 

where t denotes the specific volume, u the velocity, and a the viscosity/capillarity 
ratio, while q > and f> > are parameters. We assume here the following 
convex/concave/convex pressure law: 

P"(t)>0, t e (0,a)U(c,+oo) 
P"(t)<0, xe(fl,c), p'(a)>0. 

To tackle the traveling wave analysis we consider for definiteness a 2-wave 
issuing from (to, Uq) at — °o with speed A > 0, so 

A (t - To) + u - u = 0, 

A (u - « ) - p(t) + p(T ) = -ajS(T)|T'|V + t". 

We perform a phase plane analysis in the plane (t, t'): the equations consist of 
a second-order differential equation plus an algebraic equation. Fix a left-hand 
state To and a speed A within the interval where there exist three other equilibria 
1\, tz, ^3- These data are constrained to satisfy the following entropy inequality 

-A U' + F = -a P(t)\t'Hu') 2 < 

with entropy 
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and entropy flux 

F := up(T) + A (t') 2 + iii"-Kfl jS(t) |tT m'. 

Lemma 5.5 (Classification of equilibria). Consider the Van der Waals-type model 
above. For all q > 0, the equilibria (to, 0) and (t2, 0) are saddle points (two real 
eigenvalues with opposite signs). For q = and i = 1, 3, tte po/«f (t„ 0) z's 

• a stable node (two negative eigenvalues) ijp'(%i) + A 2 < («A/3(t,)) 2 /4, and 

• a sfaWe spira/ (two eigenvalues with the same negative real part and with opposite 
sign and non-zero imaginary parts) z/p'(T ; ) + A 2 > (aAjS(T,)) 2 /4. 

For q > Othe equilibria (ti, 0) and (T3, 0) are centers (two purely imaginary eigenval- 
ues). 

In 1 10 1, it is proven that there exists a decreasing sequence of diffusion/dispersion 
ratio a„(xo, A) — > for n > such that: 

• For a = a n there exists a nonclassical traveling wave with n oscillations 
connecting To to T2. 

• For a e {oc2 m +2/ a 2m+i) U («0/+°°)/ there exists a classical traveling wave 
connecting To to T\ . 

• For a 6 (a2m+i/«2m) there exists a classical traveling wave connecting To 

to T 3 . 

In comparison, in the case of a single inflection point one has a single critical 
value ocq(tq, A), only. Here, we have infinitely many non-monotone, nonclassical 
trajectories associated with a sequence a n — > 0. Fhis new feature observed with 
the van der Waals model indicates that the right-hand side across a given 
undercompressive wave is not unique, and several kinetic functions should be 
introduced, leading also to non-uniqueness for the Riemann problem. 

6 Existence and uniqueness theory for nonclassical 
entropy solutions 

6.1 Daf ermos' front tracking scheme 

For simplicity, consider a conservation law ll4.ll with concave-convex flux and 
impose the initial data 

u(x,Q) = u (x) u eBV(W), (6.1) 

where Uo rX is a bounded measure and the total variation TV(uq) represents the 
total mass of this measure. 
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By considering the nonclassical Rieman solver based on some kinetic func- 
tion cp" and following Dafermos [51], we construct a piecewise constant approx- 
imation u : K + X R — » ]R, as follows. First, one approximates the initial data Uq 
with a piecewise constant function w (0, •). At the time t = 0, one then solves a 
Riemann problem at each jump point of u h (0, •). If necessary, one replaces rar- 
efaction waves by several small fronts, traveling with the Rankine-Hugoniot 
speed. Then, at each interaction of waves, solve a new Riemann problem and 
continue the procedure inductively in order to construct a globally defined, 
piecewise constant approximate solution u h = u h {t,x). 

Clearly, in order to prove the convergence of the above approximation 
method, several difficulties must be overcome. First of all, one needs to show 
that the total number of wave fronts as well as the total number of interaction 
points remain finite for all fixed time. In the scalar case under consideration 
and for concave-convex flux-functions, this is actually an easy matter since each 
nonclassical Riemann solution contains at most two outgoing waves. Most 
importantly, one needs to derive a uniform bound (independent of h) on the 
total variation TV(u h (t, •)). However, due to the lack of monotonicity of the 
nonclassical Riemann solver, the (standard) total variation may increase at 
interactions. For systems, further difficulties arise due to the lack of regularity 
of the wave curves, and one must also control nonlinear interactions between 
waves of different characteristic families. 

Our assumptions on the kinetic function are very mild. We require that 
(p° : R — » K is Lipschitz continuous, monotone decreasing, and that the second 
iterate of gr is a strict contraction: for K e (0, 1) 

\<p° o (p\u)\ < K\u\, u±0. (6.2) 

Since \qft o (p^w)) < \u\ for u + 0, this is equivalent to imposing Lip !(=0 ((p 1 ' o cp b ) < 
1, that is, only a condition on nonclassical shocks with infinitesimally small 
strength. For systems of equations, similar conditions make sense and are 
realistic for the application. 

6.2 Generalized wave strength 

Following Baiti, LeFloch, and Piccoli [6] and Laforest and LeFloch l93l , we 
define the following generalized wave strength 

o{u-,u + ) := \ip(u.) - \p(u + )\, ip(u) :=\ U ' " > °' (6.3) 

[<pl(u), u<0. 

This definition has several advantages. First, it compares states with the same 
sign. Second, it is "equivalent" to the standard definition of strength, in the 
sense that _ 

C\u- - u+\ < o(u-,u+) < C\u- - u+\. 
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Third, it enjoys a continuity property as U+ crosses <p"(u_), during a transition 
from a single crossing shock to a two wave pattern: 

cr(w_, <p"(«_)) = |u_ -<Pq ° ¥*"( M -)| 

= | u - -<Po° <P b ("-)| + |<Po ° <P b ("-) - <Po ° <P*( M -)| 
= a( M _,<pV)) + ff(<pW<pV))- 

We then define a generalized total variation functional, for a piece wise 
constant function u = u(t, •) made of shock or rarefaction fronts (u", it"), by 

V(u(tj) 

a 

which again is "equivalent" to the standard total variation 

TV(u(t)) := \u a _ - u%\. 

a 

A classification of all possible wave interaction patterns is given in [3, 102], 
and about 20 cases must be distinguished. Certain cases give rise to an increase 
of the standard total variation, which is not even proportional to the (smallest) 
strength of the incoming waves. For instance, a (decreasing) classical shock 
may interact with a (decreasing) rarefaction coming from the righ-hand side, 
and transform into a (decreasing) nonclassical shock followed by an (increas- 
ing!) classical shock. Here, the outgoing wave profile is non-monotone and the 
standard total variation TV(u h (t)j increases. However, v(u h (t)^ does decrease. 
Another possible interaction is provided by a (decreasing) nonclassical shock 
which hits an (increasing) classical shock and transforms itself into a (decreas- 
ing) classical shock. In that case, significant decay of the total variation occurs, 
and both the standard total variation TV(u h (tfj and the generalized one v(u h (tfj 
decrease. 



6.3 Existence theory 

We have the following estimates and existence theory. 

Proposition 6.1 (Diminishing generalized total variation property). Consider 
a kinetic function satisfying the strict contraction property d6.2b . Then, along a se- 
quence of front tracking approximations based on the nonclassical Riemann solver, the 
generalized total variation functional V = V\u h (t)J is non-increasing. 

Theorem 6.2 (Existence of nonclassical entropy solutions (9J |93]). Consider a 
kinetic function compatible with a convex entropy U and satisfying the strict con- 
traction property j6.2i . For each initial data uq e BV(R), the wave front tracking 
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approximations u" = u h (t,x) constructed from the nonclassical Riemann solver satisfy 



l|M ft (f)llL»(R) < II«oIIl»(r), 
TV(u\t)) < TV(u ), 
\\u h (t)-u h (s)\\ LHK) <\t-s\ 



and converge in 



L 1 to a weak solution 



u = u 



(t, x) e Lip([0, +oo), L^R)) n L°°([0, +00), BV(R)) 



of the initial value problem, which satisfies the entropy inequality 



U(u) t + F(u) x < 0. 



The uniquenesss of nonclassical entropy solutions is established in Baiti, 
LeFloch, and Piccoli |8] within the class of functions with tame variation. 

We note that pre-compactness of the sequence of approximate solutions 
follows from Helly's compactness theorem. The behavior near u = is im- 
portant to prevent a blow-up of the total variation. In fact, the condition 
(p^ (0-) <p b '(0+) < 1 is indeed satisfied by kinetic functions generated by the 
nonlinear diffusion-dispersion model 



provided p < 1/2. Counter-example of blow-up of the total variation are 
available if cp b (0) = -1. See El- 
Some additional existence results are available. Perturbations of a given 
nonclassical wave are analyzed in LeFloch 11001 , Corli and Sable-Tougeron 
HZllia, Colombo and Corli EHHHl, Hattori ED, Laforest and LeFloch 11941 . 
For a version of Glimm's wave interaction potential adapted to nonclassical so- 
lutions, we refer to [93 [. The L 1 continuous dependence of nonclassical entropy 
solutions is still an open problem, and it would be interesting to generalize to 
nonclassical shocks the techniques developed by Bressan et al. [27|, LeFloch et 
al. 1Z9J M BQ21, and Liu and Yang (m) . 

Remark 6.3. An alternative strategy to establish Theorem \6J\ is developed in [7 1 which, 
however, applies only to scalar equations and requires stronger conditions on the kinetic 
function. The proposed technique of proof therein is a decomposition of the real line into 
intervals where the approximate solution is alternatively increasing/decreasing. It re- 
quires a suitable application ofFillipov-Dafermos's theory of generalized characteristics 
to track the maxima and minima of the approximate solution. 
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7 Finite difference schemes with controled dissipa- 
tion 



7.1 The role of the equivalent equation 

For simplicity in the presentation, we consider the case of the scalar conserva- 
tion law 

u t + f(u) x = e u xx + ae 2 u xxx , (7.1) 

and formulate the following question. Denoting by u a the limit when e — > 
and by qr a the associated kinetic function, can we design a numerical scheme 
converging precisely to the function u a ? 

One immediate and positive answer is provided by Glimm-type schemes, 
for which theoretical convergence results have been described in earlier sec- 
tions. This convergence was illustrated by numerical experiments performed 
with the Glimm scheme in [38 J. Another successful strategy is based on the 
level set technique implemented in [77, 44, 126] for a nonlinear elasticity model, 
with trilinear law, in two spatial dimensions, exhibiting complex interfaces with 
needles attached to the boundary. In addition, methods combining differences 
and interface tracking were also developed [156 34, 24 J, which ensure that the 
interface is sharp and (almost) exactly propagated. 

In the present section, we focus on finite difference schemes and follow 
Hayes and LeFloch (72). Let u^ x be some numerical solution and v a := 
lim^ji^o w^* be its limit. We are assuming that, at least at the practical level, 
the scheme is converging in a strong sense and does generate nonclassical 
shocks whose dynamics can be described by a kinetic function i/^. 

It was observed in [72J that 

if * u a , i/£ * (pi 

even if the scheme is conservative, consistent, high-order accurate, etc. The 
point is that in both the continuous model and the discrete scheme, small scale 
features are critical to the selection of shocks. The balance between diffusive 
and dispersive features determines which shocks are selected. Note in pass- 
ing that in nonconservative systems, the competition takes place between the 
(hyperbolic) propagation part and the (viscous) regularization [76]. 

These small scale features cannot be quite the same at the continuous and 
discrete levels, since a continuous dynamical system of ordinary differential 
equations cannot be exactly represented by a discrete dynamical system of finite 
difference equations. Consequently, finite difference schemes do not converge 
to the correct weak solution when small-scaled are the driving factor for the 
selection of shock waves. 

It was proposed [72J that i/^ should be an accurate approximation of (p\ and 
schemes with controled dissipation were developed in |[72lll091l36ll37H105lll07l , 
which rely on high-order accurate, discrete hyperbolic flux and high-order 
discretizations of the augmented terms (diffusion, dispersion). More precisely, 
it is here required that the equivalent equation of the scheme coincide with the 
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augmented physical model, up to a sufficiently high order of accuracy For 
instance, for ( |7.H , we require that after Taylor expanding the coefficients of the 
numerical scheme 

u t + f(u) x = Ax u xx + a (Ax) 2 Uxxx + 0{Axf, (7.2) 

for p > 3 at least. 

It was conjectured by the author that as p — > 00 the kinetic function i/^ p 
associated with a scheme having the equivalent equation i7.2) converges to the 
exact kinetic function that is, 

Support for this conjecture was recently provided by LeFloch and Mohamadian 
11071 , who performed extensive numerical tests for several models including 
the generalized Camassa-Holm and Van der Waals ones. 

7.2 The role of entropy conservative schemes 

In addition, the role of entropy conservative schemes was stressed in 172111091 . 
Schemes have been built from higher order accurate, entropy conservative, 
discrete flux. Such schemes satisfy discrete versions of the physically relevant 
entropy inequality, hence preserve exactly (and globally in time) an approxi- 
mate entropy balance 

The design of entropy conservative schemes is based on entropy variable. 
Consider a system of conservation laws CQ} endowed with an entropy pair 
(U,F). Suppose that U strictly convex or, more generally, f(u) can be expressed 
as a function of v, and let v(u) = VU(m) 6 V := V!i(U) be the entropy variable. 
Finally, set f(u) = g(v), F(u) = G(v), and B(v) = Dg(v). It is easily checked that 
B(v) is symmetric, since Dg(v) = D/(m)D 2 !J(m)" 1 . So, there exists \p(v) such that 
g = Vi/> and, in fact, 

ip(v) = v ■ g{v) - G(v). 

On a regular mesh Xj = jh (j = ...,-1,0,1,...), consider (2p + l)-point, 
conservative, semi-discrete schemes 

dt U ' ~~ h [g i +1/2 g hW'' 
where Uj = uAt) represents an approximation of u(xj, t). The discrete flux 

S/+1/2 = g*( v i-p+i> "•> v i+v)> v i = V!J (";) 
must be consistent with the exact flux g, i.e. 

g\v,...,v) = g(v). 

Tadmor in [146 1 introduced the notion of entropy conservative schemes 
and focused on essentially three-point schemes, for which g*{v- v +\, ■ ■ ■ ,v v ) = 
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g(v) when Vq = v \ = v, which implies second-order accuracy. The two-point 
numerical flux [146[ 

g*{v Q ,v x )= I g(v + s (v x - v )) ds, v ,vieU. 
Jo 

yields an entropy conservative scheme, satisfying 

jU{u,) + \{G) +ll2 -G)_ V2 ) = 0, 

with 

1 1 
G*(v ,vx) =- (G(v ) + G(vx)) + -(v + Vx)g'(v ,Vx) 

1 

This scheme admits the following second-order accurate, (conservative) equiv- 
alent equation 111091 

u t + f(u) x = — ( - g(v) xx + -v x - Dg(v) x ) x 

with v — Vlf(w). This is not sufficient for our purpose of tackling the diffusion- 
dispersion model. 

High-order entropy conservative schemes were discovered by LeFloch and 
Rohde 111091 . Later, generalizations to arbitrarily high order were found in 
I105H1471H481 . 

Theorem 7.1 (Third-order entropy conservative schemes). Consider the semi- 
disrete (2p + V)-point scheme 

Jt U i = ~l ( £; +1 /2 " S'j-xnl gj+i/z = g*(Pj- P +i> ■■■ ■ V i+p) 

with numerical flux g* defined as follows, from any symmetric N x N matrices 
B*(v- p+2 , ■■■ , v p ), 

g*(v- v+ i,--- ,v p ) 

= j" g(vo +s(vi- v fj ds - i( (v 2 - Vi) ■ B\v- v+2 , ■ ■ • , v p ) 

- (v - v-x) ■ B*(v- p+ x, ■■■ , Vp-x)). 
This scheme is entropy conservative for the entropy If, with entropy flux 
G*{v- v +x,--- ,v p ) 

1 1 

= 2^0 + Vx) • g\V-p+l, ■ ■ ■ , Vp) - -{lp*(V- p+ 2, ■ ■ ■ , Vp) + 4>*{V- V+ X, ■■■ , Vp-x)). 
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and 

il>*(v- p +2, ■ • ■ , v p ) 

1 

= vi ■ g(vt) - G{vi) + — (v x - v ) ■ B'(v- p+2/ ■■■ , Vp)(vt - v 2 ). 

When p = 2 and B*{v, v, v) = B(v)[ = Dg(vfj, this five-point scheme is third-order, at 
least. 

Although schemes with controled dissipation based on an analysis of their 
equivalent equation do not converge to the exact solution determined by a 
given augmented model, still they provide a large class of practically useful 
schemes and allow one to ensure that the numerical kinetic function approaches 
the exact one. 

It is now established numerically that kinetic functions exist and are mono- 
tone for a large class of physically relevant models including thin liquid films, 
generalized Camassa-Holm, nonlinear phase transitions, van der Waals fluids 
(for small shocks), and magnetohydrodynamics. 

Computing the kinetic function has been found to be very useful to in- 
vestigate the effects of the diffusion/dispersion ratio, regularization, order of 
accuracy of the schemes, the efficiency of the schemes, as well as to make 
comparisons between several physical models. 

In addition, it should be noted that kinetic functions may be associated with 
schemes 1721 : The Beam- Warming scheme (for concave-convex flux) produces 
non-classical shocks, while no such shocks are observed with the Lax-Wendroff 
scheme. All of this depends crucially on the sign of the numerical dispersion 
coefficient. 

As mentioned earlier, similar issues arise in dealing with the numerical 
approximation of nonlinear hyperbolic systems in nonconservative form for 
which we refer to [76, 33] and the references therein. 

8 Concluding remarks 

Let us conclude by mentioning a few more issues of interest. 

Vanishing diffusion-dispersion limits in the context of the initial value prob- 
lem have been investigated by several approaches. Tartar's compensated com- 
pactness method was applied to treat one-dimensional scalar conservation 
laws, by Schonbek [138[, Hayes and LeFloch [71 J, and LeFloch and Natalini 
|108[. The 2x2 system of nonlinear elasticity system was tackled in 1751 . 
Singular limits for the Camassa-Holm equation were analyzed by Coclite and 
Karlsen l40l . In all these works, conditions are imposed on the diffusion and 
dispersion parameters which are mild enough to allow for nonclassical shocks 
in the limit. 

Another strategy based on DiPerna's measure-valued solutions [56] applies 
to multidimensional conservation laws. Strong convergence results were estab- 
lished by Correia and LeFloch [49, 50] and Kondo and LeFloch [90]. The gen- 
eralization to discontinuous flux is provided in Holden, Karlsen, and Mitrovic 
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175] . DiPerna's theorem requires all of the entropy inequalities, and therefore 
these results do not cover the regime of parameters allowing for nonclassical 
shocks. 

A third approach is based on Lions, Perthame, and Tadmor's kinetic for- 
mulation 1 118 1 and again applies to multidimensional conservation laws. It 
was first observed by Hwang and Tzavaras [82] that the kinetic formulation 
extends to singular limits to conservation laws when diffusive and dispersive 
parameters are kept in balance. See also Hwang [80 . 81] and Kwon |92|. The 
kinetic formulation was recently extended to non-local regularizations |89|. 

In another direction, LeFloch and Shearer |110| introduced a nucleation 
criterion and a Riemann solver with kinetic and nucleation for scalar conserva- 
tion laws. They could cope with problems in which the traveling wave analysis 
does not select a unique solution to the Riemann problem when, even after 
imposing a kinetic relation for undercompressive shocks, one is still left with 
a classical and a nonclassical Riemann solution. Earlier on, Abeyaratne and 
Knowles [ 1 1 had introduced a nucleation criterion in the context of the model 
of elastodynamics with a trilinear equation of state. 

Roughly speaking, the nucleation criterion imposes that a "sufficiently 
large" initial jump always "nucleates". The qualitative properties of the Rie- 
mann solver with kinetic and nucleation were investigated in 11101 : prescribing 
the set of admissible waves does not uniquely determine the Riemann solution, 
and instability phenomena such as "splitting-merging" wave structures take 
place. The analysis was recently extended to hyperbolic systems of conserva- 
tion laws [94J and was also further investigated numerically [115] . It is expected 
that the nucleation criterion will be particulary relevant for higher-order regu- 
larizations such as the one in the thin liquid film model. 

In conclusion, a large class of scalar equations and 2x2 systems are now 
well-understood. The Riemann problem is uniquely determined by an entropy 
inequality and a kinetic relation, and the kinetic function can be determined 
by an analysis of traveling solutions. For these example, the Riemann solution 
depends continuously upon its initial data. However, for other models such 
as the thin film model or the hyperbolic-elliptic model of phase transitions, the 
existing theory has limitations, and this suggests challenging open problems. 

Due to a lack of space (and time), still many other issues could not be 
treated in this review. Kinetic relations are relevant also for nonconservative 
hyperbolic systems [19], and should play an important role in the mathematical 
modeling of multi-fluid and turbulence models. For an extensive literature on 
the nonlinear stability of undercompressive shock waves (including for multi- 
dimensional problems), we refer to works by T.-P. Liu, Metivier, Williams, 
Zumbrun, and others. See [73] [H [122] |T23] Q33) [137) and the references 
therein. 



36 



Acknowledgments 



I am particularly grateful to Constantine M. Dafermos, Tai-Ping Liu, and 
Michael Shearer for scientific discussions on this research work and for their 
encouragments to pursue it. 

This paper was written as part of the international research program on 
Nonlinear Partial Differential Equations held at the Centre for Advanced Study 
of the Norwegian Academy of Science and Letters during the Academic Year 
2008-09. I am very grateful to Helge Holden and Kenneth Karlsen for their invi- 
tation and hospitality and for the opportunity to give a short-course, originally 
entitled "Small-scale dependent shock waves. Theory, approximation, and 
applications". The author was also supported by a DFG-CNRS collaborative 
grant between France and Germany on "Micro-Macro Modeling and Simula- 
tion of Liquid- Vapor Flows", as well as by the Centre National de la Recherche 
Scientifique (CNRS) and the Agence Nationale de la Recherche (ANR) via the 
grant 06-2-134423. 

References 

[1] R. Abeyaratne and J.K. Rnowles, Kinetic relations and the propagation of phase 
boundaries in solids, Arch. Rational Mech. Anal. 114 (1991), 119-154. 

[2] R. Abeyaratne and J.K. Knowles, Implications of viscosity and strain-gradient 
effects for the kinetics of propagating phase boundaries in solids, SIAM J. Appl. 
Math. 51 (1991), 1205-1221. 

[3] D. Amadori, P. Baiti, PG. LeFloch and B. Piccoli, Nonclassical shocks and the 
Cauchy problem for nonconvex conservation laws, J. Differential Equations 151 
(1999), 345-372. 

[4] P. Antonelli and M. Marcati, On the finite energy weak solutions to a system in 
quantum fluid dynamics, Comm. Math. Phys. 287 (2009), 657-686. 

[5] A. Azevedo, D. Marchesin, B.J. Plohr, and K. Zumbrun, Non-uniqueness of solu- 
tions of Riemann problems caused by 2-cycles of shock waves, Proc. Fifth Internat. 
Conf. on Hyperbolic Problems: theory, numerics, applications, J. Glimm, M.J. 
Graham, J.W. Grove, and B.J. Plohr, ed., World Scientific Editions, 1996, pp. 43-51. 

[6] P. Baiti, PG LeFloch, and B. Piccoli, Nonclassical shocks and the Cauchy problem. 
General conservation laws, Contemporary Math. 238 (1999), 1-25. 

[7] P. Baiti, PG LeFloch, and B. Piccoli, BV Stability via generalized characteristics for 
nonclassical solutions of conservation laws, EQUADIFF'99, Proc. Internat. Conf. 
Differ. Equ., Berlin, August 1999, B. Fiedler, K. Groger, and J.Sprekels, editors, 
World Sc. Publ., River Edge, NY, 2000, pp. 289-294. 

[8] P. Baiti, PG. LeFloch, and B. Piccoli, Uniqueness of classical and nonclassical 
solutions for nonlinear hyperbolic systems, J. Differential Equations 172 (2001), 
59-82. 



37 



[9] P. Baiti, P.G. LeFloch, and B. Piccoli, Existence theory for nonclassical entropy 
solutions: scalar conservation laws, Z. Angew. Math. Phys. 55 (2004), 927-945. 

[10] N. Bedjaoui, C. Chalons, F. Coquel, and P.G. LeFloch, Non-monotonic traveling 
waves in van der Waals fluids, Anal. Appl. 3 (2005), 419^46. 

[11] N. Bedjaoui and P.G. LeFloch, Diffusive-dispersive traveling waves and kinetic 
relations I. Non-convex hyperbolic conservation laws, J. Differential Equations 178 
(2002), 574-607. 

[12] N. Bedjaoui and P.G. LeFloch, Diffusive-dispersive traveling waves and kinetic 
relations II. A hyperbolic-elliptic model of phase transitions dynamics, Proc. Royal 
Soc. Edinburgh 132A (2002), 545-565. 

[13] N. Bedjaoui and P.G. LeFloch, Diffusive-dispersive traveling waves and kinetic 
relations III. An hyperbolic model from nonlinear elastodynamics, Ann. Univ. 
Ferrara Sc. Mat. 47 (2002), 117-144. 

[14] N. Bedjaoui and P.G. LeFloch, Diffusive-dispersive traveling waves and kinetic 
relations. IV. Compressible Euler system, Chinese Ann. Appl. Math. 24 (2003), 
17-34. 

[15] N. Bedjaoui and P.G. LeFloch, Diffusive-dispersive traveling waves and kinetic 
relations. V. Singular diffusion and dispersion terms, Proc. Royal Soc. Edinburgh 
134A (2004), 815-844. 

[16] S. Benzoni-Gavage, Stability of subsonic planar phase boundaries in a van der 
Waals fluid, Arch. Ration. Mech. Anal. 150 (1999), 23-55. 

[17] E. Beretta, J. Hulshof, and LA. Peletier, On an ordinary differential equation from 
forced coating flow, J. Differential Equations 130 (1996), 247-265. 

[18] C. Berthon and F. Coquel, Nonlinear projection methods for multi-entropies 
Navier-Stokes systems, Math, of Comp. 76 (2007), 1163-1194. 

[19] C. Berthon, F. Coquel, and P.G. LeFloch, Kinetic relations for nonconservative 
hyperbolic systems and applications, unpublished notes (2002). 

[20] A. Bertozzi, A. Munch, and M. Shearer, Undercompressive shocks in thin film flow, 
Phys. D 134 (1999), 431-464. 

[21] A. Bertozzi, A. Munch, M. Shearer, and K. Zumbrun, Stability of compressive and 
undercompressive thin film traveling waves, European J. Appl. Math. 12 (2001), 
253-291. 

[22] A. Bertozzi and M. Shearer, Existence of undercompressive traveling waves in thin 
film equations, SIAM J. Math. Anal. 32 (2000), 194-213. 

[23] T. Bohme, W. Dreyer, and W.H. Miiller, Determination of stiffness and higher 
gradient coefficients by means of the embedded-atom method, Contin. Mech. 
Thermodyn. 18 (2007), 411-441. 

[24] B. Boutin, C. Chalons, F. Lagoutiere, and P.G. LeFloch, Convergent and conserva- 
tive schemes for nonclassical solutions based on kinetic relations, Interfaces and 
Free Boundaries 10 (2008), 399^21. 



38 



[25] S. Bianchini and A. Bressan, Vanishing viscosity solutions of nonlinear hyperbolic 
systems, Ann. of Math. 161 (2005), 223-342. 



[26] F. Bouchut, C. Klingenberg, and K. Waagan, A multiwave approximate Riemann 
solver for ideal MHD based on relaxation. I. Theoretical framework, Numer. Math. 
108 (2007), 7-42. 

[27] A. Bressan, Hyperbolic systems of conservation laws. The one-dimensional Cauchy prob- 
lem, Oxford Lecture Series Math. Appl., 20, Oxford University Press, Oxford, 2000. 

[28] A. Bressan and A. Constantin, Global dissipative solutions of the Camassa-Holm 
equation, Anal. Appl. 5 (2007), 1-27. 

[29] A. Bressan and P.G. LeFloch, Uniqueness of entropy solutions for systems of con- 
servation laws, Arch. Rational Mech. Anal. 140 (1999), 301-331. 

[30] M. Brio and J.K. Hunter, Rotationally invariant hyperbolic waves, Comm. Pure 
Appl. Math. 43 (1990), 1037-1053. 

[31] G. Caginalp and X. Chen, Convergence of the phase field model to its sharp 
interface limits, European J. Appl. Math. 9 (1998), 417-445. 

[32] G. Carboux and B. Hanouzet, Relaxation approximation of some initial-boundary 
value problem for p-systems, Commun. Math. Sci. 5 (2007), 187-203. 

[33] M.J. Castro, P.G. LeFloch, M.L. Munoz-Ruiz, and C. Pares, Why many theories of 
shock waves are necessary. Convergence error in formally path-consistent schemes, 
J. Comput. Phys. 227 (2008), 8107-8129. 

[34] C. Chalons, Transport-equilibrium schemes for computing nonclassical shocks. 
Scalar conservation laws, Numer. Methods Partial Differential Equations 24 (2008), 
1127-1147. 

[35] C. Chalons and F. Coquel, Numerical capture of shock solutions of nonconser- 
vative hyperbolic systems via kinetic functions. Analysis and simulation of fluid 
dynamics, Adv. Math. Fluid Mech., pp. 45-68, Birkhauser, Basel, 2007. 

[36] C. Chalons and P.G. LeFloch, High-order entropy conservative schemes and kinetic 
relations for van der Waals fluids, J. Comput. Phys. 167 (2001), 1-23. 

[37] C. Chalons and P.G. LeFloch, A fully discrete scheme for diffusive-dispersive 
conservation laws, Numerische Math. 89 (2001), 493-509. 

[38] C. Chalons and P.G. LeFloch, Computing undercompressive waves with the ran- 
dom choice scheme : nonclassical shock waves, Interfaces and Free Boundaries 5 
(2003), 129-158. 

[39] M. Charlotte and L. Truskinovsky, Towards multi-scale continuum elasticity theory, 
Contin. Mech. Thermodyn. 20 (2008), 133-161. 

[40] G.M. Coclite and K.H. Karlsen, A singular limit problem for conservation laws 
related to the Camassa-Holm shallow water equation, Comm. Partial Differential 
Equations 31 (2006), 1253-1272. 



39 



[41] R.M. Colombo and A. Corli, Continuous dependence in conservation laws with 
phase transitions, SIAM J. Math. Anal. 31 (1999), 34-62. 

[42] R.M. Colombo and A. Corli, Stability of the Riemann semigroup with respect to 
the kinetic condition, Quart. Appl. Math. 62 (2004), 541-551. 

[43] R.M. Colombo and A. Corli, Sonic and kinetic phase transitions with applications 
to Chapman-Jouguet deflagration, Math. Methods Appl. Sci. 27 (2004), 843-864. 

[44] F. Coquel, D. Diehl, C. Merkle, and C. Rohde, Sharp and diffuse interface methods 
for phase transition problems in liquid-vapour flows, in "Numerical methods for 
hyperbolic and kinetic problems", pp. 239-270, IRMA Lect. Math. Theor. Phys., 7, 
Eur. Math. Soc, Zurich, 2005. 

[45] A. Corli and M. Sable-Tougeron, Kinetic stabilization of a nonlinear sonic phase 
boundary, Arch. Rational Mech. Anal. 152 (2000), 1-63. 

[46] A. Corli and M. Sable-Tougeron, Stability of contact discontinuities under pertur- 
bations of bounded variation, Rend. Sem. Mat. Univ. Padova 97 (1997), 35-60. 

[47] A. Corli and M. Sable-Tougeron, Perturbations of bounded variation of a strong 
shock wave, J. Differential Equations 138 (1997), 195-228. 

[48] A. Corli and M. Sable-Tougeron, Kinetic stabilization of a nonlinear sonic phase 
boundary, Arch. Ration. Mech. Anal. 152 (2000), 1-63. 

[49] J.M. Correia and PC LeFloch, Nonlinear diffusive-dispersive limits for multidi- 
mensional conservation laws, in "Advances in Nonlinear P.D.E.'s and Related Ar- 
eas", (Beijing, 1997), World Sci. Publ., River Edge, NJ, 1998, pp. 103-123. Available 
at: http://arxiv.org/abs/0810.1880 

[50] J.M. Correia and PC LeFloch, Nonlinear hyperbolic conservation laws, in "Non- 
linear evolution equations and their applications" (Macau, 1998), World Sci. Publ., 
River Edge, NJ, 1999, pp. 21^4. 

[51] CM. Dafermos, Polygonal approximations of solutions of the initial value problem 
for a conservation law, J. Math. Anal. Appl. 38 (1972), 33^1. 

[52] CM. Dafermos, Hyperbolic systems of conservation laws, Proceedings "Systems 
of Nonlinear Partial Differential Equations", J.M. Ball editor, NATO Adv. Sci. Series 
C, 111, Dordrecht D. Reidel (1983), 25-70. 

[53] CM. Dafermos, Hyperbolic conservation laws in continuum physics, Grundlehren 
Math. Wissenschaften Series 325, Springer Verlag, 2000. 

[54] G. Dal Maso, PC LeFloch, and F. Murat, Definition and weak stability of noncon- 
servative products, J. Math. Pures Appl. 74 (1995), 483-548. 

[55] C. De Lellis, F. Otto, and M. Westdickenberg, Minimal entropy conditions for 
Burgers equation, Quart. Appl. Math. 62 (2004), 687-700. 

[56] R.J. DiPerna, Measure-valued solutions to conservation laws, Arch. Rational Mech. 
Anal. 88 (1985), 223-270. 



40 



[57] W. Dreyer and M. Herrmann, Numerical experiments on the modulation theory 
for the nonlinear atomic chain, Phys. D 237 (2008), 255-282. 

[58] W. Dreyer and M. Kunik, Maximum entropy principle revisited, Contin. Mech. 
Thermodyn. 10 (1998), 331-347. 

[59] W. E and D. Li, The Andersen thermostat in molecular dynamics, Comm. Pure 
Appl. Math. 61 (2008), 96-136. 

[60] H.T. Fan and H.-L. Liu, Pattern formation, wave propagation and stability in 
conservation laws with slow diffusion and fast reaction, J. Differ. Hyper. Equa. 1 
(2004), 605-636. 

[61] H.-T. Fan and M. Slemrod, The Riemann problem for systems of conservation 
laws of mixed type, in "Shock induced transitions and phase structures in general 
media", Workshop held in Minneapolis (USA), Oct. 1990, Dunn J.E. (ed.) et al., 
IMA Vol. Math. Appl. 52 (1993), pp. 61-91. 

[62] H. Freistuhler, Dynamical stability and vanishing viscosity: A case study of a 
non-strictly hyperbolic system of conservation laws, Comm. Pure Appl. Math. 45 
(1992), 561-582. 

[63] H. Freistuhler and E.B. Pitman, A numerical study of a rotationally degenerate 
hyperbolic system. I. The Riemann problem, J. Comput. Phys. 100 (1992), n306- 
321. 

[64] H. Freistuhler and E.B. Pitman, A numerical study of a rotationally degenerate 
hyperbolic system. II. The Cauchy problem, SIAM J. Numer. Anal. 32 (1995), 741- 
753. 

[65] H. Frid and I-S. Liu, Phase transitions and oscillation waves in an elastic bar, 
in"Fourth Workshop on Partial Differential Equations", Part I (Rio de Janeiro, 
1995), Mat. Contemp. 10 (1996), 123-135. 

[66] P. Goatin and PC LeFloch, Sharp L 1 continuous dependence of solutions of 
bounded variation for hyperbolic systems of conservation laws, Arch. Ration. 
Mech. Anal. 157 (2001), 35-73. 

[67] M. Grinfeld, Nonisothermal dynamic phase transitions, Quart. Appl. Math. 47 
(1989), 71-84. 

[68] H. Hattori, The Riemann problem for a van der Waals fluid with entropy rate 
admissibility criterion — isothermal case, Arch. Rational Mech. Anal. 92 (1986), 
247-263. 

[69] H. Hattori, The existence and large time behavior of solutions to a system related 
to a phase transition problem, SIAM J. Math. Anal. 34 (2003), 774-804. 

[70] H. Hattori, Existence of solutions with moving phase boundaries in thermoelas- 
ticity, J. Hyperbolic Differ. Equ. 5 (2008), 589-611 

[71] B.T Hayes and PC. LeFloch, Nonclassical shocks and kinetic relations. Scalar 
conservation laws, Arch. Rational Mech. Anal. 139 (1997), 1-56. 



41 



[72] B.T. Hayes and P.G. LeFloch, Nonclassical shocks and kinetic relations. Finite 
difference schemes, SIAM J. Numer. Anal. 35 (1998), 2169-2194. 

[73] B.T. Hayes and P.G. LeFloch, Nonclassical shock waves and kinetic relations. 
Strictly hyperbolic systems, SIAM J. Math. Anal. 31 (2000), 941-991. (Preprint # 
357, CMAP, Ecole Polytechnique, Palaiseau, France, Nov. 1996.) 

[74] B.T. Hayes and M. Shearer, Undercompressive shocks and Riemann problems for 
scalar conservation laws with non-convex fluxes, Proc. Roy. Soc. Edinburgh Sect. 
A 129 (1999), 733-754. 

[75] H. Holden, K.H. Karlsen, and D. Mitrovic, Zero diffusion-dispersion smoothing 
limits for a scalar conservation law with discontinuous flux function, International 
Journal of Differential Equations (2009). 

[76] T.Y. Hou and P.G. LeFloch, Why nonconservative schemes converge to wrong 
solutions. Error analysis, Math, of Comput. 62 (1994), 497-530. 

[77] T.Y. Hou, P. Rosakis, and P.G. LeFloch, A level set approach to the computation of 
twinning and phase transition dynamics, J. Comput. Phys. 150 (1999), 302-331. 

[78] P. Howard and K. Zumbrun, Stability of undercompressive shock profiles, J. Dif- 
ferential Equations 225 (2006), 308-360. 

[79] J. Hu and P.G. LeFloch, L 1 continuous dependence property for systems of conser- 
vation laws, Arch. Ration. Mech. Anal. 151 (2000), 45-93. 

[80] S. Hwang, Kinetic decomposition for the generalized BBM-Burgers equations with 
dissipative term, Proc. Roy. Soc. Edinburgh Sect. A 134 (2004), 1149-1162. 

[81] S. Hwang, Singular limit problem of the Camassa-Holm type equation, J. Differ- 
ential Equations 235 (2007), 74-84. 

[82] S. Hwang and A. Tzavaras, Kinetic decomposition of approximate solutions to con- 
servation laws: application to relaxation and diffusion-dispersion approximations, 
Comm. Partial Differential Equations 27 (2002), 1229-1254. 

[83] T. Iguchi and P.G. LeFloch, Existence theory for hyperbolic systems of conservation 
laws with general flux-functions, Arch. Rational Mech. Anal. 168 (2003), 165-244. 

[84] E. Isaacson, D. Marchesin, C.F. Palmeira, and B.J. Plohr, A global formalism for 
nonlinear waves in conservation laws, Comm. Math. Phys. 146 (1992), 505-552. 

[85] D. Jacobs, W.R. McKinney, and M. Shearer, Traveling wave solutions of the modi- 
fied Korteweg-deVries Burgers equation, J. Differential Equations 116 (1995), 448- 
467. 

[86] J.W. Jerome, The mathematical study and approximation of semi-conductor mod- 
els. Large-scale matrix problems and the numerical solution of partial differential 
equations (Lancaster, 1992), pp. 157-204, Adv. Numer. Anal., Ill, Oxford Univ. 
Press, New York, 1994. 

[87] K.T Joseph and P.G. LeFloch, Singular limits in phase dynamics with physical 
viscosity and capillarity, Proc. Royal Soc. Edinburgh 137A (2007), 1287-1312. 



42 



[88] K.T. Joseph and P.G. LeFloch, Singular limits for the Riemann problem: general 
diffusion, relaxation, and boundary conditions, C.R. Math. Acad. Sci. Paris 344 
(2007), 59-64. 

[89] F. Kissling, P.G. LeFloch, and C. Rohde, Singular limits and kinetic decomposition 
for a non-local diffusion-dispersion problem, J. Differential Equations 247 (2009), 
3338-3356. 

[90] C. Kondo and P.G. LeFloch, Zero diffusion-dispersion limits for hyperbolic con- 
servation laws, SIAM Math. Anal. 33 (2002), 1320-1329. 

[91] S.N. Kruzkov, First order quasilinear equations with several independent variables, 
Mat. Sb. (NS.) 81 (123) (1970), 228-25. 

[92] Y.-S. Kwon, Diffusion-dispersion limits for multidimensional scalar conservation 
laws with source terms, J. Differential Equations 246 (2009), 1883-1893. 

[93] M. Laforest and P.G. LeFloch, Diminishing functionals for nonclassical entropy 
solutions selected by kinetic relations, Port. Math. (2010). 

[94] M. Laforest and P.G. LeFloch, in preparation. 

[95] P.D. Lax, Hyperbolic systems of conservation laws, II, Comm. Pure Appl. Math. 
10 (1957), 537-566. 

[96] P.D. Lax, Hyperbolic systems of conservation laws and the mathematical theory of shock 
waves, Regional Confer. Series in Appl. Math. 11, SIAM, Philadelphia, 1973. 

[97] P.D. Lax and CD. Levermore, The small dispersion limit of the Korteweg-deVries 
equation, Comm. Pure Appl. Math. 36 (1983), 253-290. 

[98] P.G. LeFloch, Entropy weak solutions to nonlinear hyperbolic systems in noncon- 
servative form, Comm. Part. Diff. Equa. 13 (1988), 669-727. 

[99] P.G. LeFloch, An existence and uniqueness result for two nonstrictly hyperbolic 
systems, in "Nonlinear evolution equations that change type", pp. 126-138, IMA 
Vol. Math. Appl, 27, Springer, New York, 1990. 

[100] P.G. LeFloch, Propagating phase boundaries: formulation of the problem and 
existence via the Glimm scheme, Arch. Rational Mech. Anal. 123 (1993), 153-197. 

[101] P.G. LeFloch, An introduction to nonclassical shocks of systems of conservation 
laws, International School on Hyperbolic Problems, Freiburg, Germany, Oct. 97, 
D. Kroner, M. Ohlberger and C. Rohde eds., Lect. Notes Comput. Eng., Vol. 5, 
Springer Verlag, 1999, pp. 28-72. 

[102] P.G. LeFloch, Hyperbolic systems of conservation laws: the theory of classical and 
nonclassical shock waves, Lecture Notes in Mathematics, E.T.H. Zurich, Birkhauser, 
2002. 

[103] P.G. LeFloch, Graph solutions of nonlinear hyperbolic systems, J. Hyperbolic 
Differ. Equ. 1 (2004), 643-689. 

[104] P.G. LeFloch and T.-P. Liu, Existence theory for nonlinear hyperbolic systems in 
nonconservative form, Forum Math. 5 (1993), 261-280. 



43 



[105] P.G. LeFloch, J.-M. Mercier, and C. Rohde, Fully discrete entropy conservative 
schemes of arbitrary order, SIAM J. Numer. Anal. 40 (2002), 1968-1992. 

[106] P.G. LeFloch and S. Mishra, in preparation. 

[107] P.G. LeFloch and M. Mohamadian, Why many shock wave theories are necessary. 
Fourth-order models, kinetic functions, and equivalent equations, J. Comput. Phys. 
227 (2008), 4162-4189. 

[108] P.G. LeFloch and R. Natalini, Conservation laws with vanishing nonlinear diffu- 
sion and dispersion, Nonlinear Analysis 36 (1999), 213-230. 

[109] P.G. LeFloch and C. Rohde, High-order schemes, entropy inequalities, and non- 
classical shocks, SIAM J. Numer. Anal. 37 (2000), 2023-2060. 

[110] P.G. LeFloch and M. Shearer, Nonclassical Riemann solvers with nucleation, Proc. 
Royal Soc. Edinburgh 134A (2004), 941-964. 

[Ill] P.G. LeFloch and M.D. Thanh, Nonclassical Riemann solvers and kinetic relations. 
III. A nonconvex hyperbolic model for Van der Waals fluids, Electr. J. Diff. Equa. 
72 (2000), 1-19. 

[112] P.G. LeFloch and M.D. Thanh, Nonclassical Riemann solvers and kinetic relations 

I. A nonconvex hyperbolic model of phase transitions, Z. Angew. Math. Phys. 52 
(2001), 597-619. 

[113] P.G. LeFloch and M.D. Thanh, Nonclassical Riemann solvers and kinetic relations. 

II. An hyperbolic-elliptic model of phase transitions, Proc. Royal Soc. Edinburgh 
132A (2001), 181-219. 

[114] P.G. LeFloch and M.D. Thanh, Properties of Rankine-Hugoniot curves for Van der 
Waals fluids, Japan J. Indust. Applied Math. 20 (2003), 211-238. 

[115] R. Levy and M. Shearer, Comparison of two dynamic contact line models for 
driven thin liquid films, European J. Appl. Math. 15 (2004), 625-642. 

[116] R. Levy and M. Shearer, Kinetics and nucleation for driven thin film flow, Phys. 
D 209 (2005), 145-163. 

[117] H.-L. Li and R.-H. Pan, Zero relaxation limit for piecewise smooth solutions to a 
rate-type viscoelastic system in the presence of shocks, J. Math. Anal. Appl. 252 
(2000), 298-324. 

[118] P.-L. Lions, B. Perthame, and E. Tadmor, A kinetic formulation of multidimen- 
sional scalar conservation laws and related equations, J. Amer. Math. Soc. 7 (1994), 
169-191. 

[119] T.-P. Liu, The Riemann problem for general 2x2 conservation laws, Trans. Amer. 
Math. Soc. 199 (1974), 89-112. 

[120] T.-P. Liu, Admissible solutions of hyperbolic conservation laws, Mem. Amer. Math. Soc. 
30, 1981. 

[121] T.-P. Liu and T. Yang, Well-posedness theory for hyperbolic conservation laws, 
Comm. Pure Appl. Math. 52 (1999), 1553-1586. 



44 



[122] T.-P. Liu and K. Zumbrun, Nonlinear stability of an undercompressive shock for 
complex Burgers equation, Comm. Math. Phys. 168 (1995), 163-186. 

[123] T.-P. Liu and K. Zumbrun, Nonlinear stability of general undercompressive shock 
waves, Comm. Math. Phys. 174 (1995), 319-345. 

[124] J.-M. Mercier and B. Piccoli, Global continuous Riemann solver for nonlinear 
elasticity, Arch. Ration. Mech. Anal. 156 (2001), 89-119. 

[125] J.-M. Mercier and B. Piccoli, Admissible Riemann solvers for genuinely nonlinear 
p-systems of mixed type, J. Differential Equations 180 (2002), 395-426. 

[126] C. Merkle and C. Rohde, The sharp-interface approach for fluids with phase 
change: Riemann problems and ghost fluid techniques, M2AN Math. Model. 
Numer. Anal. 41 (2007), 1089-1123. 

[127] A. Munch, Shock transitions in Marangoni gravity-driven thin-film flow, Nonlin- 
earity 13 (2000), 731-746. 

[128] S.-C. Ngan and L. Truskinovsky L., Thermo-elastic aspects of dynamic nucleation, 
J. Mech. Phys. Solids 50 (2002), 1193-1229. 

[129] O. Oleinik, Discontinuous solutions of nonlinear differential equations, Transl. 
Amer. Math. Soc. 26 (1963), 95-172. 

[130] F. Otto and M. Westdickenberg, Convergence of thin film approximation for a 
scalar conservation law, J. Hyperbolic Differ. Equ. 2 (2005), 183-199. 

[131] E.Y. Panov, On the theory of entropy solutions of the Cauchy problem for a class of 
nonstrictly hyperbolic systems of conservation laws, Sb. Math. 191 (2000), 121-150. 

[132] E.Y. Panov, On the theory of generalized entropy sub- and supersolutions of the 
Cauchy problem for a first-order quasilinear equation, Differ. Equ. 37 (2001), 272- 
280. 

[133] M. Raoofi and K. Zumbrun, Stability of undercompressive viscous shock profiles 
of hyperbolic-parabolic systems, J. Differential Equations 246 (2009), 1539-1567. 

[134] A. Ratz and A. Voigt, A diffuse-interface approximation for surface diffusion 
including adatoms, Nonlinearity 20 (2007), 177-192. 

[135] C. Rohde, Scalar conservation laws with mixed local and non-local diffusion- 
dispersion terms, SIAM J. Math. Anal. 37 (2005), 103-129. 

[136] C. Rohde, On local and non-local Navier-Stokes-Korteweg systems for liquid- 
vapour phase transitions, ZAMM Z. Angew. Math. Mech. 85 (2005), 839-857. 

[137] S. Schecter and M. Shearer, Undercompressive shocks for nonstrictly hyperbolic 
conservation laws, Dynamics Diff. Equa. 3 (1991), 199-271. 

[138] M.E. Schonbek, Convergence of solutions to nonlinear dispersive equations, 
Comm. Part. Diff. Equa. 7 (1982), 959-1000. 

[139] S. Schulze and M. Shearer, Undercompressive shocks for a system of hyperbolic 
conservation laws with cubic nonlinearity, J. Math. Anal. Appl. 229 (1999), 344-362. 



45 



[140] M. Shearer, The Riemann problem for the planar motion of an elastic string, J. 
Differential Equations 61 (1986), 149-163. 

[141] M. Shearer, Dynamic phase transitions in a van der Waals gas, Quart. Appl. Math. 
46 (1988), 631-636. 

[142] M. Shearer and Y. Yang, The Riemann problem for a system of conservation laws 
of mixed type with a cubic nonlinearity, Proc. Roy. Soc. Edinburgh Sect. A 125 
(1995), 675-699. 

[143] M. Slemrod, Admissibility criteria for propagating phase boundaries in a van der 
Waals fluid, Arch. Rational Mech. Anal. 81 (1983), 301-315. 

[144] M. Slemrod, A limiting viscosity approach to the Riemann problem for materials 
exhibiting change of phase, Arch. Rational Mech. Anal. 105 (1989), 327-365. 

[145] I. Suliciu, On modelling phase transition by means of rate-type constitutive equa- 
tions, shock wave structure, Int. J. Ing. Sci. 28 (1990), 827-841. 

[146] E. Tadmor, The numerical viscosity of entropy stable schemes for systems of 
conservation laws. I, Math. Comp. 49 (1987), 91-103. 

[147] E. Tadmor, Entropy stability theory for difference approximations of nonlinear 
conservation laws and related time-dependent problems, Acta Numer. 12 (2003), 
451-512. 

[148] E. Tadmor and W. Zhong, Entropy stable approximations of Navier-Stokes equa- 
tions with no artificial numerical viscosity, J. Hyperbolic Differ. Equ. 3 (2006), 
529-559. 

[149] L. Truskinovsky, Dynamics of non-equilibrium phase boundaries in a heat con- 
ducting nonlinear elastic medium, J. Appl. Math, and Mech. (PMM) 51 (1987), 
777-784. 

[150] L. Truskinovsky, Kinks versus shocks, in "Shock induced transitions and phase 
structures in general media", R. Fosdick, E. Dunn, and M. Slemrod ed., IMA Vol. 
Math. Appl., Vol. 52, Springer- Verlag, New York (1993), pp. 185-229. 

[151] L. Truskinovsky, Transition to detonation in dynamic phase changes, Arch. Ratio- 
nal Mech. Anal. 125 (1994), 375-397. 

[152] L. Truskinovsky and A. Vainchtein, Quasicontinuum models of dynamic phase 
transitions, Contin. Mech. Thermodyn. 18 (2006), 1-21. 

[153] C.J. Van Duijn, L.A. Peletier, and I.S. Pop, A new class of entropy solutions of the 
Buckley-Leverett equation, SIAM J. Math. Anal. 39 (2007), 507-536. 

[154] A.I. Volpert, The space BV and quasi-linear equations, Mat. USSR Sb. 2 (1967), 
225-267. 

[155] B. Wendroff The Riemann problem for materials with non-convex equations of 
state. I. Isentropic flow, J. Math. Anal. Appl. 38 (1972), 454-466. 

[156] X.-G. Zhong, TY Hou, and PG. LeFloch, Computational methods for propagating 
phase boundaries, J. Comput. Phys. 124 (1996), 192-216. 



46 



